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ABSTRACT 


Over the years, there has been significant variation in the filmwise steam 
condensation data at NPS on horizontal low-integral finned tubes. With a view to 
increasing the accuracy of the data, inserts were used inside the tubes to reduce 
inside thermal resistance; however, significant discrepancies then occurred in the 
calculated outside coefficient when compared to data taken without an insert. These 
discrepancies arose due to the data reduction technique which assumes a known 
inside heat-transfer resistance and subtracts this from a measured overall resistance. 
If the assumed value on the inside is inaccurate, then the outside value is equally 
inaccurate. 

The present work uses an instrumented smooth tube to obtain accurate inside 
heat-transfer correlations both with and without inserts and uses these to obtain 
accurate outside coefficients for a family of uninstrumented finned tubes with a view 


to finding an optimum fin spacing for steam condensation. 
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I. INTRODUCTION 


A. BACKGROUND 

A reduction in size and weight of all types of heat exchangers aboard Naval 
vessels will allow more efficient use of space. The benefits might include greater 
equipment accessibility for maintenance or greater heat exchanger capacity (without 
a corresponding increase in size and weight) with a corresponding increase in fuel 
efficiency. 

For the past ten years, the Naval Postgraduate School in collaboration with the 
David Taylor Research Center and the National Science Foundation has conducted 
research that is directed at the development of smaller, more efficient steam 
condensers. Improved designs can result in significant space savings, always a 
primary concern on Naval vessels, especially submarines. 

Uncertainties in past data using steam were apparently due to the lack of 
detailed information about the inside heat-transfer correlations used during the data 
reduction process. Previously, the standard Sieder-Tate correlation was assumed to 
be valid for the inside heat-transfer coefficient, but it may not be the best correlation 
to use with the particular test arrangement used in this research program. 

A large amount of enhanced condensation data has been collected in previous 
studies at the Naval Postgraduate School on more than 90 different condenser tubes 


of varying fin height, fin spacing, and tube material, with a view to finding an 


optimum fin geometry for both steam and refrigerant condensation. Satisfactory 
results have been obtained with refrigerant data using R-113. However, some 
troublesome questions of possible contamination and instrument inaccuracy still 
remain with the steam data. It is felt that in order to address these questions, a 
fundamental evaluation of the heat transfer apparatus on which this data was 
collected and of the data reduction process was considered appropriate. Current 
data, taken on a carefully cleaned and calibrated apparatus, could be compared to 
previously recorded data and a determination as to its validity and reproducibility 
could be made. In addition, a thorough evaluation of the best inside heat-transfer 


correlation would lead to more reliable steam condensation results. 


B. CONDENSATION 

Condensation occurs when a vapor is cooled below its saturation temperature, 
or when a vapor/gas mixture is cooled below its dew point. Surface condensation 
occurs in condensers when a cooled surface (kept at a temperature below the 
saturation temperature of the vapor) contacts the vapor. The vapor molecules that 
contact such a surface stick to that surface and condense into liquid molecules. 
Condensation may occur in one of the following modes: filmwise, dropwise, or mixed 
mode (a combination of filmwise and dropwise) condensation. In the filmwise mode, 
the liquid wets the cold surface to form a continuous film. If the liquid does not wet 


the surface but instead forms discrete drops on the cold surface, dropwise or mixed 


mode condensation will occur and is often caused by some form of contamination. 
[Ref. 1] 

The condensate forming on the tube surface offers a resistance to heat transfer 
between the vapor and the surface, which increases with the thickness of the liquid 
layer. Even though dropwise condensation results in much larger heat-transfer 
coefficients than filmwise condensation, it is difficult to maintain a stable dropwise 
condition over prolonged periods. Therefore, in most cases condenser design 
calculations are based on the assumption of filmwise condensation, resulting in lower 
heat-transfer coefficients and more conservative designs. [Ref. 1] 

In a condenser, the coolant side, tube wall, and vapor side thermal resistances 
control the heat transfer rate from vapor to coolant. Also, for experimental work 
we always use clean tubes, but in real condensers tubeside fouling can play an 
important role an increasing the coolant side resistance. The magnitudes of these 
resistances depend on the fluid, tube geometry, and flow conditions on the vapor and 
coolant side. For steam condensation, it is the coolant side thermal resistance which 
tends to dominate. Methods to lower this inside resistance include the use of inserts 
or roped tubing to promote turbulence, thereby raising the convection heat-transfer 
coefficient. However, such modifications lead to increased pressure drop through 
the tubes, which must be compensated for by providing extra pumping capacity. 
Heat transfer through the tube wall is conductive and is fixed once tube thickness 


and material are selected. The vapor side resistance is due to the condensate film 


which forms on the outside of the tube. For filmwise condensation, the outside 


resistance can be reduced by the addition of low integral fins. These have the effect 
of not only increasing the outside surface area of the tube, but also of thinning the 
condensate film around the fins due to surface tension forces. Too small a fin 
spacing may result in condensate flooding, whereas too large a fin spacing 
approaches the smooth tube case; there should be an optimum fin spacing 
somewhere in between these two extremes. Horizontal fin spacing is therefore of 
prime importance, and finding the optimum spacing is one of the objectives of this 


long-term research program. 


C. CONDENSATION RESEARCH AT NAVAL POSTGRADUATE SCHOOL 

The research effort at NPS has included the study of differing fin dimensions 
(i.e. fin height, fin width, fin spacing) on low-integral finned horizontal tubes. 
Experimentation has included the use of three different test fluids (steam, R-113, 
and ethylene glycol) under various operational conditions using a number of different 
tube diameters. 

Van Petten [Ref. 2] provides a summary of research efforts at NPS through the 
end of 1988. Van Petten and subsequent researchers have analyzed small, medium, 
and large diameter finned tubes to find the optimum fin spacing for maximum heat- 
transfer enhancement of the fluids mentioned above. However, discrepancies found 
by Guttendorf [Ref. 3] in the data processing technique (modified Wilson plot), 


which resulted in different values of heat-transfer enhancement (for the same tube 





under the same operating conditions depending on whether an insert was or was not 
used), have raised doubts about the accuracy of the inside heat-transfer correlation. 

Rouk [Ref. 4] investigated the use of an optimization technique to predict the 
inside heat-transfer correlation. When the optimization effort proved unsuccessful, 
he next used the instrumented smooth tube data of Georgiadis [Ref. 5] to develop 
an inside heat transfer correlation, but could not find a correlation with sufficient 
accuracy based on previous data. He recommended that once an overhaul on the 
test apparatus was complete, an increase in data precision would allow the 
development of an accurate inside heat-transfer correlation. This work is a follow 
on effort to develop inside heat-transfer correlations which can predict the value of 
the inside heat-transfer coefficient with good accuracy under a variety of flow 
conditions. Once an inside correlation can be found, the object of this effort is to 
reprocess previous data and see if the discrepancies reported by former researchers 


for finned tubes on this apparatus can be rectified. 


D. OBJECTIVES 
The main objectives of this thesis were to: 
1. Disassemble and meticulously clean the apparatus to eliminate any existing 
contamination with a view to eliminating dropwise condensation problems 


experienced in the past. 


2. Carefully reassemble the apparatus using new gasket material, and make 
modifications to improve system performance. 


3. Recalibrate all system instrumentation to ensure the greatest achievable 
accuracy. 


Investigate the possibility of manufacturing large, medium, and small diameter 
instrumented smooth tubes. 


Use the new instrumented tubes and the one existing medium diameter 
smooth instrumented tube (fabricated by Poole [Ref. 6]) to obtain accurate 
inside heat-transfer correlations for a number of insert types as well as the 
no insert condition. 


Evaluate the accuracy of the currently used data processing technique 
(modified Wilson plot) using instrumented tube data. 


Reprocess previous data using the new correlations with a view to comparing 
current and past smooth tube and finned tube data to provide continuity with 
previous studies. 


П. LITERATURE SURVEY 


A. INTRODUCTION 

When a vapor condenses in the filmwise mode on a smooth horizontal tube it 
forms a thin continuous film of condensate on the surface of the tube. The 
condensate film thickens around the tube due to gravity. This condensate film 
provides a resistance to heat transfer which may be lowered through the use of fins. 
For quite some time, it was thought to be impractical to use finned tubes with high 
surface tension fluids such as water, due to condensate retention and flooding 
between the fins. However, a number of studies conducted on finned tubes using 
steam have shown that substantial heat-transfer enhancement may be achieved. 

A significant amount of research at the Naval Postgraduate School and 
elsewhere has addressed the issue of optimum fin height, thickness, and spacing 
required for maximum heat transfer. Yau et al [Ref. 7] reported that "with an 
increase in fin density, up to a limit (this limit is not yet known in a generalized 
manner), the heat-transfer coefficient increases at a rate faster than the increase in 
the outside area due to the presence of fins. This additional enhancement is due to 
the thinning effect of the surface-tension forces on the condensate film. 
Unfortunately, surface-tension forces also adversely affect heat transfer by causing 


condensate to be retained between fins" [Ref. 8]. Katz et al [Ref. 9] also found that 


on finned tubes the portion of the surface occupied by condensate is dependent upon 
the ratio of condensate surface-tension to density and the fin geometry. 
Condensate retention and the behavior of the condensate film on the tube 
surface under various conditions are critical parameters in the heat transfer process 
on horizontal finned tubes. Several models have been developed to predict this 
behavior and the reader is referred to an extensive review of horizontal finned tube 


heat transfer by Marto [Ref. 10] for a more detailed coverage of the topic. 


B. VAPOR SIDE CONSIDERATIONS 

The filmwise condensation of vapor on a horizontal tube is a complex two- 
phase heat transfer process, for which a suitably complex model would be required 
to accurately predict heat transfer performance under all conditions. 

In 1916, Nusselt [Ref. 11] set forth his theoretical work on the study of laminar 
filmwise condensation of a "stationary" vapor on a vertical or inclined plate and a 


horizontal tube. Nusselt’s simplifying assumptions included the following [Ref. 12]: 


1. Pure saturated vapor 

2. Negligible vapor velocity (U,=0) 

3. Heat transfer across the condensate film by conduction only 

4. Laminar condensate flow governed only by gravitational and viscous forces 
5. Condensate properties constant | 

6. Isothermal condensing surface 


7. Negligible interface temperature drop 


Nusselt's result for the mean heat-transfer coefficient for a horizontal tube was 














obtained: 
3 1/4 
ми =0.728 7 РАР PVR Me (2.1) 
ром) 
or 
1/3 
kp kp, 
Nu =0.655| ҰРАР PPR (2.2) 
u,D,q 

where: 

Nu = mean Nusselt number 

К, = thermal conductivity of condensate film (W/m - k) 

р; = condensate film density (kg/m’) 

р, = vapor density (kg/m’) 

g = gravitational constant (9.81 m/s’) 

h, - specific enthalpy of vaporization (J/kg) 

HM, = dynamic viscosity of condensate film (N-s/m*) 

D,  - outside tube diameter (m) 

AT; = average temperature difference across condensate film (K) 

q = heat flux based on outside area (Q/A,) (W/m?) 


Many workers have improved on Nusselt’s theoretical analysis by accounting 
for some of the terms he neglected through his simplifying assumptions. However, 
equations (2.1) and (2.2) have been found to be remarkably accurate over a wide 
range of conditions for a stationary vapor. High vapor velocity can increase film 
condensation heat transfer substantially. This enhancement, which refers to the 
amount of heat transfer above or below the value predicted by the Nusselt analysis, 


is due to the effect of thinning the condensate film. However, vapor shear is the one 


assumption which if applied can lead to significant increases in the heat-transfer 
coefficient. 
The theoretical result of Shekriladze and Gomelauri (1966) [Ref. 13], who 


considered interfacial shear stress due to vapor velocity, is shown in equation (2.3). 





0641 411.9 (2.3) 
20 
where: 
Nu = Nusselt number for the vapor side 
Ве. = two phase Reynolds number, (p, U, Р/и,) 


For steam condensation, the empirically derived correlation of Fujii et al [Ref. 


14] is shown in equation (2.4). 


Nu 


Ке, 





= 0.96F 5 (2.4) 


The Nusselt expression (equation (2.1) can be expressed in similar form: 


Nu 


Ке, 





=0.728F 4 (25) 


Whereas the vapor velocity, U,, cancels out in the Nusselt expression 
(stationary vapor assumption), the Fujii correlation includes the vapor velocity effect. 


Therefore we can expect equation (2.4) to more accurately predict steam-side heat 


10 


transfer coefficients for those cases where vapor velocity begins to have a significant 
impact. 
For further review of basic theoretical studies on the subject of laminar film 


condensation on smooth tubes the reader is referred to Rose [Ref. 15]. 


C. COOLANT SIDE CONSIDERATIONS 
For a turbulent flow regime inside a pipe (Re > 10,000), a number of coolant- 
side correlations have been used; many of these have taken the form: 
Nu - C,Re" Pr" (2.6) 


where: 


Nu mean coolant Nusselt number for turbulent flow 
C; = correlating coefficient 

Re coolant Reynolds number 

Pr coolant Prandtl number 


The most common correlations with the same form as equation (2.6) are that 


of Dittus and Boelter (1930) (Ref. 16]: 


Nu =0.023 Re®® pr°4 (2.7) 


and Colburn (1933) [Ref. 17]: 


Nu - 0.023 Re 9? pr !'? (2.8) 


11 


Sieder and Tate (1936) [Ref. 18] applied a correction factor to equation (2.8) 
to account for cases in which the bulk to inner wall temperature difference is large 


enough to cause substantial variations in coolant viscosity as follows: 


Nu - 0.027 Re95 pr (y, jy. 914 (2.10) 
where: 
H. = coolant viscosity evaluated at mean bulk temperature (N:s/m?) 
ш, = coolant viscosity evaluated at mean inner tube wall temperature 


(N:s/m^) 


The fluid properties in equations (2.7), (2.8), and (2.9) are evaluated at mean 


coolant bulk temperature T, bux 


Tu UOTA (2.11) 
where: 
T, = tube coolant inlet temperature (K) 
T, = tube coolant outlet temperature (K) 


Equations (2.7), (2.8) and (2.9) are valid for Re» 10! and 0.7 < Рг < 100, апа 
were developed for long smooth pipes with no inserts [Ref. 12]. 

The use of inserts and the effect of bends close to the tube entrance region 
can affect the values of both the leading coefficient, C; , and Reynolds number 
exponent, m, in equation (2.6). One of the major focuses of this study is to 
determine the values of C; and m for the no insert, wire wrap insert, and Heatex 


insert cases. 
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Other well-known turbulent pipe flow heat-transfer correlations (i.e. Petukhov- 
Popov [Ref. 19], Sleicher-Rouse [Ref. 20], etc.) and the results of an ANL (Argonne 
National Laboratory) [Ref. 21] study which evaluated several such correlations for 


accuracy are reviewed in section VI C. 
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ШІ. APPARATUS AND SYSTEM INSTRUMENTATION 


A. SYSTEM OVERVIEW 

The apparatus used for this research was basically the same as was used by Van 
Petten [Ref. 2] and Guttendorf [Ref. 3] with certain modifications. A system 
schematic is provided in Figure 1. Steam generated from distilled water in the .30 
m diameter pyrex glass boiler using ten 4 kW, 440 V Watlow immersion heaters was 
the working medium for this set of experiments. From the boiler section the steam 
passed up through a reducing section and a 2.13 m straight length of pyrex glass 
piping, (ID of 0.15 m), it was then turned through 180 degrees using two 90 degree 
pyrex glass elbows, and then descended down a 1.52 m straight length of pyrex glass 
piping. Тһе steam then entered the stainless steel test section containing the 
horizontally mounted condenser tube (see Figures 1 and 2); any steam not 
condensing there was condensed in the auxiliary condenser located just beneath the 
test section. The auxiliary condenser was constructed of a single copper coil mounted 
to a stainless steel baseplate enclosed within a pyrex glass condenser section. 
Coolant flow through the auxiliary condenser was used to control system pressure, 
and all condensed liquid was returned via the condenser baseplate drain to the boiler 
section by gravity. 

Coolant for the auxiliary condenser was provided via a throttled water 


connection with associated flowmeter. Coolant flow through the single horizontal 
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Figure 1. Schematic of Single Tube Test Apparatus 


15 


13714 
“ЕР 


“om 7 
МОЛАМ 2 
all A: омн5л 3319 Др 832945 
— = — — ER ETT 
| ANS 2 
u 
2 








тізе 
3" 3gOr 381 


Figure 2. Schematic of Test Section Insert 
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tube was provided by a separate system consisting of a sump tank with two 
centrifugal pumps connected in series. Coolant flow rate was measured by a 
carefully calibrated flowmeter. By varying coolant flow-rate through the single 
horizontal tube, the rate of steam condensation on the tube (and hence heat-transfer 
coefficient) could be varied. 

Non-condensible gases were removed using the vacuum pump system shown in 
Figure 3. The condensing coil for this purge system, located in the sump tank, 
served to condense steam carried through the vacuum line during the purging 
process. The vacuum line took its suction from the base of the auxiliary condenser, 
the coolest spot in the apparatus and the place where non-condensible gases (1.e. air) 


were most likely to accumulate. 


B. SYSTEM INSTRUMENTATION 

The power to the 440 V heaters was controlled through a panel mounted 
potentiometer. A description of the power calculation for input into the data 
acquisition system can be found in Poole [Ref. 6]. 


System pressure was monitored in three ways: 


1. A Setra model 204 pressure transducer 
2. A Heise solid front pressure gauge (visual reading only) 


3. System saturation pressure from vapor temperature measurement 
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Figure 3. Schematic of Purging System and Cooling Water Sump 
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System vapor temperature was monitored using both teflon, and metal sheathed 
type-T copper/constantan thermocouples located juxtaposed in the test section; this 
position was just upstream of the test condenser tube. Condensate temperature was 
also monitored using a teflon coated type-T copper/constantan thermocouple located 
on the condensate return line between the auxiliary condenser and boiler. Coolant 


temperature rise in the condenser tube was measured using four methods: 


1. Two teflon coated type-T copper/constantan thermocouples 
2. Two metal sheathed type-T copper/constantan thermocouples 
3. Two Hewlett-Packard 2804A quartz crystal thermometers 


4. A ten-junction teflon coated type-T copper/constantan thermopile 


These were all placed at the inlet to and exit from the condenser tube; at the 
outlet, all thermocouples were placed just downstream of a coolant mixing chamber. 

Two data reduction programs were used to collect and reduce data on this 
apparatus; DRPINST', and "DRPKS". The instrumented tube constructed by Poole 
[Ref. 6] was used to determine an accurate inside heat-transfer correlation for inserts 
used; this instrumented tube contained six wall thermocouples. For the instrumented 
tube the appropriate calibration equations were accessed in the data acquisition 
program "DRPINST". For non-instrumented tubes, the data reduction program 
"DRPKS" was used. Fluid property 2. used in the data reduction programs 


are given in Appendix A and calibrations were conducted for all system 
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instrumentation (flowmeter, thermocouples, pressure transducer, etc.) and are 
included in Appendix B. 

]he data as monitored by the aforementioned system instrumentation was 
processed by an HP-3497A data acquisition system controlled by an HP-9826A 
computer provided with the correct data acquisition program. The raw data was 
processed and stored on computer disks. Program channel assignments are given in 


Table 1. 


C. SYSTEM MODIFICATIONS 
At the beginning of this investigation the apparatus was entirely disassembled 
to facilitate complete overhaul of the system. Modifications and details of assembly 


were as follows: 


1. The apparatus was taken apart piece by piece, inspected, meticulously cleaned 
with a warm solution of Sparkleen biodegradable soap and subjected to a 
complete acetone rinse prior to reassembly. 


2. A new pyrex glass riser section above the boiler was 0.31 m shorter than the 
previous section and allowed the addition of a new aluminum stand on which 
to place the heater baseplate. This new stand allowed much easier access to 
the 440 V heater wiring plus the adjustable legs allowed level adjustment of 
the entire apparatus, to ensure proper alignment of the single horizontal tube. 


3. The two pyrex glass elbows were replaced. 
4. Every gasket in the system was replaced (using Buna-N rubber) and, using a 


standard star torque pattern, all flanged joints were tightened to a final 
torqued of 30 inch-pounds (manufacturer recommended maximum torque 
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Table 1. DATA ACQUISITION SYSTEM CHANNEL ASSIGNMENT 


CHANNEL 


40 


4] 


42 


43 
44 
45 


46 
47 
48 


49 


50 


51 


52 


53 


60 


61 
62 
64 


QUANTITY 
MEASURED 
Vapor Temperature 
Vapor Temperature 


Room Temperature 


Tube Coolant in 1 
Tube Coolant out 1 
Tube Coolant in 2 


Tube Coolant out 2 
Condensate Return 


Instrumented Tube 
Instrumented Tube 
Instrumented Tube 
Instrumented Tube 
Instrumented Tube 
Instrumented Tube 
10 - Junction 
Thermopile 
Voltage 


Current 


Pressure Transducer 
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APPLICABLE TYPE 
DESIGNATION 


T-57 
(See Appendix B.1) 


1-56 
(See Appendix B.1) 


T-58 
(See Appendix B.1) 


1-58 
1-58 
T-55 

(See Appendix B.1) 
T-55 
Т-58 


Т-57 ог Т-55 
as applicable 


1-57 ог Т-55 
as applicable 


1-57 or T-55 
as applicable 


1-57 or 1-55 
as applicable 


Т-57 ог Т-55 
as applicable 


1-57 or 1-55 
as applicable 


Т-61 
(See Appendix B.1) 


NA 
МА 
(See Appendix B.4) 


specification was 60 inch-pounds). The bolts holding the flanges together 
should be checked for tightness on a regular basis, since the thermal cycling 
of the apparatus has been shown to cause loosening of nut and bolt 
assemblies. 


The previous vacuum pump was replaced with a Gast model 2567-V108 
vacuum pump which could draw vacuum to 130 mmHg with an installed check 
valve to prevent pump back spin when the pump was stopped. The new 
pump could draw vacuum much more rapidly, but could not operate at less 
than 130 mmHg. Remaining non-condensible gases were removed by flushing 
the system with steam. The steam flushing procedure for removal of non- 
condensible gases is given in the operating procedures section. Once the 
system is completely filled with steam, operating pressures well below 130 
mmHg could be achieved utilizing the auxiliary condenser. 


The double coil auxiliary condenser was replaced with a single coil used 
originally by Van Petten [Ref. 2]. The single coil was not coated with the 
special oxide coating used by Guttendorf [Ref. 3]. It was felt to be superior 
to the double coil in that the baseplate welds were of much higher quality and 
were preferred on the basis of vacuum tightness. 


The aluminum side plates attached to the pyrex glass auxiliary condenser 
housing were replaced with new stainless steel side plates with penetrations 
for pressure bleed, vacuum line, and a pressure transducer. These three 
penetrations were fitted with screw threaded stainless steel connectors. The 
stainless steel connectors were heli-arc welded in place. Prior to this 
modification (completed 24 January 1992) a leak test conducted from 21 
December 199] to 2 January 1992 revealed a mean vacuum leak rate of ~ 3.4 
mmHg per day (see Figure C.1). A subsequent leak test conducted from 6 
February 1992 to 19 February 1992 showed an improvement in the mean leak 
rate to —1.7 mmHg per day (see Figure C2). 


System instrument modifications included the addition of the Setra pressure 
transducer and the Heise pressure gauge, and the removal of the mercury 
manometer. All system instrumentation was recalibrated and the results 
incorporated into the data reduction programs. 


Finally the apparatus was lagged with Halstead insulating foam to reduce heat 


loss as much as possible. The test section, which was left uncovered 
previously, was also lagged. 
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А. 


IV EQUIPMENT OPERATION AND EXPERIMENTAL PROCEDURE 


SYSTEM STARTUP AND SHUTDOWN PROCEDURES 
Startup of the system is accomplished in the following manner: 


Ensure distilled water level in the boiler is 4 to 6 inches above the top of the 
heating elements. The boiler is filled by gravity drain via a hose connection 
from the distilled water tank to the boiler fill valve. Ensure the vent valve on 
the side of the auxiliary condenser is open when filling or draining the boiler. 
The boiler may be drained by removing the hose connection and opening the 
fill valve, which allows drainage into the trench directly beneath the boiler. 


Once the boiler is filled to the appropriate level, shut the boiler fill valve and 
the distilled water tank valve. 


Shut the system vent valve. 


Turn on the data acquisition system, computer and printer. Load the 
appropriate program (DRPINST or DRPKS) and check for proper operation. 
Then check all thermocouple outputs, by stepping through the appropriate 
data acquisition system channels, to verify that all are registering ambient 
temperature. 


Open the fill valve to the coolant water sump tank to a level such that the 
tank overflow drain box does not overflow (the valve is located between the 
boiler control panel and heat pipe apparatus). 


Turn on the cooling water supply pumps and adjust the tube flow rate from 
20% to 60% of the rotameter setting and check for leaks. Reset flow rate to 
desired level. 


Open valves from tap water system to auxiliary condenser and adjust coolant 
flow rate to at least 30% and check for leaks. Reset flow rate to at least 
10%. 


Energize heaters and adjust voltage to approximately 50 volts (40 volts if the 
system is already at vacuum below 100 mmHg to limit the vibrational shock 
to the system from oversized vapor bubble formation). To energize the 
heaters there are three switches which must be placed in the on position. The 
first is located in power panel p5 located in the main hallway adjacent to the 
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10. 


11. 


12. 


lab and is labeled switch 3 / heater controller room 106. The second is the 
heater load bank circuit breaker located on the side of the boiler control 
panel. The third is the condensing rig boiler power switch on the front of the 
boiler control panel. Increase the voltage gradually in 10 volt increments to 
the desired level. 


Turn on the vacuum pump and open the vacuum line valve. Allow the 
vacuum pump to run until system pressure is below 3 psi, then shut the 
vacuum line valve just prior to turning off the vacuum pump. 


As system warmup continues and pressure increases to above 4 psi energize 
the vacuum pump as necessary to flush the non-condensible gases out of the 
system through the vacuum line by forcing the gases out with steam. То 
ensure that the non-condensibles gather at the base of the auxiliary 
condenser, where the vacuum line suction is located, ensure that the 
horizontal tube is not supplied with coolant flow, and adjust coolant flow 
through the auxiliary condenser as necessary to ensure steam is filling the 
entire system. The auxiliary condenser may be touched lightly by hand, along 
its entire length, to ensure the system is completely filled with steam; any cool 
spots indicate the presence of non-condensible gases which means that the 
flushing process is not complete. The flushing process takes 15 to 30 minutes 
to accomplish, and should be repeated periodically for long periods of 
operation. 


At the conclusion of the flushing process, shut the vacuum line valve and 
secure the vacuum pump. 


In order to ensure that filmwise condensation occurs on the tube, coolant flow 
through the tube must be initiated as follows: 
a. Allow the apparatus vapor temperature measurement (channel 40) to 


reach at least 3800 microvolts. 


b. Cut in the auxiliary coolant flow (50% or 60% level) to cool the vapor 
temperature to roughly 3200 microvolts. 


c. Secure coolant flow through the auxiliary condenser, and allow the 
vapor temperature level to climb to about 3700-3800 microvolts, which 


allows a steam blanket to cover the tube. 


d.  Initiate coolant flow through the single horizontal tube at the 80% level. 


24 


13. 


e. Cut in coolant flow to the auxiliary condenser to control pressure, and 
observe the condensation process to ensure that a condensate film has 
formed on the tube. 


Run software program DRPINST for an instrumented tube (DRPKS for 
uninstrumented tubes) by pressing "run" on the keyboard. 


To take data for an instrumented tube, the questions for DRPINST can be 


answered as follows: 


Select fluid ... Enter 0 for steam 

select option ... Enter 1 to take new data 

Enter month, date and time ... Press enter 

Enter input mode ... Enter 1 for new data 

Give a name for the raw data file ... Enter name 
Enter geometry code ... Enter 1 for finned, 0 for plain 


Select insert type ... Enter 0 for none, 1 for twisted tape, 2 for wire wrap, 3 for 
Heatex 


No. of thermocouples in wall? ... Enter 4, 5, or 6 depending on the tube 
Select tube diameter type ... Enter 2 for medium 

Enter pressure condition ... Enter 0 for vacuum, 1 for atmospheric 

Give a name for the wall temperature file ... Enter name 

Select input ... Enter 1 for short, 2 for long, or 3 for raw data 


Like to check NG (non-condensible gas) concentration ... Enter 1 for yes, 2 for 
no; you must answer yes for the first data point 


Enter flowmeter reading (%) ... Enter 2 digit number (i.e. 20 or 58 etc.) 


Connect voltage line ... Flip the voltage line toggle switch, located on the power 
contro] panel, to the on position and press enter 
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Disconnect voltage line ... Flip the voltage line toggle switch off and press enter 
Enter pressure gauge reading (Pga) ... Enter reading off gauge in psi 


Select measurement ... Enter 0 for teflon, 1 for metal sheath, 2 for quartz, 
3 for thermopile 


Change TCOOL rise? ... Enter 1 for yes, 2 for no 
OK to store this data set? ... Enter 1 for yes, 0 Ғог по 


Will there be another run? ... Enter 1 for yes, 0 for no; starts at check NG 
concentration for following runs. 


To take data for an uninstrumented tube the questions for DRPKS can be 


answered as follows: 


Select fluid ... Enter 0 for steam 

Select option ... Enter 0 to take new data 
Enter month, date and time ... Press enter 
Enter disk number ... Enter number 
Enter input mode ... Enter 0 for new data 


Select Ci ... Enter 0 to find a Ci value, 2 to use a Ci value stored in the 
program. 


Give a name for the raw data file ... Enter name 
Enter geometry code ... Enter 1 for finned, 0 for plain 


Enter insert type ... Enter 0 for none, 1 for twisted type, 2 for wire wrap, 3 for 
Heatex | 


Select tube type ... Enter 0 for thick wall (only thick wall tubes were tested) 


Select material code ... Enter 0 for copper (only copper tubes were tested) 
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15. 


Select tube diameter type ... Enter 1 for medium (no small or large diameter 
tubes were tested) 


Enter pressure condition ... Enter 0 for vacuum, 1 for atmospheric 
Want to create a file for NR vs F? ... Enter 1 for yes, 0 for no 


Give a name for plot data file ... Enter name; easiest to use the raw data file 
name preceded by a P 


Select output ... Enter 0 for short, 1 for long, 2 for raw data 


Like to check NG concentration ... Enter 1 for yes, 2 for no; you must answer 
yes for the first data point 


Enter flowmeter reading (%) ... Enter 2 digit number (i.e. 20 or 60 etc.) 
Connect voltage line ... Flip the voltage line toggle switch on and press enter 
Disconnect voltage line ... Flip the voltage line toggle switch off and press enter 
Enter pressure gauge reading (Pga) ... Enter reading off gauge in psi 


Select measurement ... Enter 0 for teflon, 1 for metal sheath, 2 for quartz, 
3 for thermopile 


Change TCOOL rise ... Enter 1 for yes, 2 for no 
OK to store this data set ... Enter 1 for yes, 0 for no 
Will there be another run ... Enter 1 for yes, 0 for no; starts at check NG 


concentration for following runs. 


Only answer the program questions up to "Enter flowmeter readings’. 
Monitor system temperature using the vapor thermocouple voltage reading 
(the program automatically resets to channel 40) closely until system warmup 
is complete. 


Monitor system temperature and pressure carefully to prevent a system over 
pressure during warmup (especially at atmospheric conditions). 
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17. 


18. 


19. 


20) 


If conducting a vacuum run, gradually adjust voltage to 90 volts (usually in 10 
volt increments). Obtain the desired operating condition by manually 
controlling coolant flow through the auxiliary condenser until channel 40 
reads 1970 + 20 microvolts (~48°C). Vapor velocity ~ 2 m/s. 


If conducting an atmospheric run, gradually adjust voltage to 175 volts from 
the 90 volt level in 10-20 volt increments. Again the desired operating 
condition is obtained by manually controlling coolant flow through the 
auxiliary condenser until channel 40 reads 4280 + 20 microvolts (~ 100°C). 
Vapor velocity ~ 1 m/s. 


Monitor the condensation process using the glass viewing window periodically 
to ensure that filmwise condensation is maintained. To clear the viewing 
window of fog and moisture increase coolant flow through the auxiliary 
condenser briefly to 50% or 60%, then reset to desired flow rate. 


When taking readings be sure to check the flowmeter setting prior to entering 
it into the computer (it has a tendency to fluctuate slightly). 


If conducting vacuum and atmospheric runs on the same day always conduct 
the vacuum run first. If the atmospheric run is done first it takes too long for 
the system to cool down to vacuum operating temperatures. 


The system is secured in the following manner: 


Secure power to the heating elements. 


Secure coolant flow through the tube, through the auxiliary condenser, and 
to the sump tank. 


If desired to maintain the system at vacuum conditions until the next run the 
shutdown is complete. Continued cooling water circulation may be used to 


assist in cooling down the system. 


To bring the system back to atmospheric conditions slowly open the vent 
valve. 


The data acquisition system may be turned off whenever it is not necessary 
to monitor system parameters. 


Periodically change distilled water in the boiler. 
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7. If an emergency should arise such as abrupt overpressurization or breakage, 
immediately secure power to the heaters and open the vent valve, then let the 
system cool down before checking the apparatus for damage. 

B. EXPERIMENTAL PROCEDURES AND OBSERVATIONS 

Water is a poor wetting medium and therefore great care was taken to ensure 
that uniform filmwise condensation was the only condensation mode occurring during 
a data run. Even though the apparatus was meticulously cleaned (as mentioned 
previously), a continuing problem with dropwise condensation manifested itself. 
Subsequent to steam cleaning the system with a Sparkleen soap solution, by 
operating the system with a soapy solution in the boiler (the solution bubbled 
through the entire apparatus), dropwise condensation was observed on the installed 
instrumented tube. After taking some data when in the dropwise condition, the tube 
was removed and rigorously cleaned using a warm Alconox soap solution with a 
scrub brush. However, after observing the filmwise mode initially, the condensation 
mechanism soon transitioned to mixed mode and then back to the dropwise mode. 

Since only a filmwise condition over several hours would suffice, the tube 
chemical treatment procedure used by Guttendorf [Ref. 3] and several other 
researchers was used to produce filmwise condensation. The tube was chemically 
treated prior to installation as follows: 

1. Clean the internal and especially de external surfaces of the tube using a soft 


brush and mild soap (using the Alconox detergent in warm water), rinse with 
acetone then rinse thoroughly with distilled water. Repeat the cleaning 
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procedure until the distilled water rinse perfectly wets the tube surface; any 
breaks in the wetting film at this point are likely to result in dropwise 
condensation spots once the tube is installed in the apparatus. 

2. Place the tube in a steam bath. 

3. Mix equal amounts of ethyl alcohol and a 50% by weight solution of sodium 
hydroxide. Keep the solution warm so that a watery consistency is 
maintained. 

4. Apply the solution to the tube with a small paint brush, retaining the tube in 
the steam bath. If the tube has not been treated previously, apply a coating 
of the solution every 10 minutes for an hour. If the tube has been previously 
treated, apply a coating every 5 minutes for a period of 20 minutes. 

5. Remove the tube from the steam bath and thoroughly rinse the tube with 
distilled water to remove any excess solution. Install the tube in the test 
section immediately, being careful not to touch the tube surface. Oil or dirt 
from any source may contaminate the tube surface and result in mixed mode 
or dropwise condensation. 

The oxide layer which forms on the tube is very thin, and has negligible thermal 
resistance and high wetting characteristics. 

Once the tube was installed in the apparatus (with the desired insert in place), 
the system startup procedure outlined in section IV A was followed to take data at 
desired conditions. 

At vacuum conditions, when single tube coolant flow was initiated with vapor 
velocity at — 2 m/s, the condensation on the tube did not develop as a perfect film 
but instead left patches where the film was broken. These patches, or streaks 
seemed to occur at regular intervals, and it was postulated that they were due to 


vortex shedding of vapor around the tube. Therefore, the procedure in the startup 


section IV A, step 12, was used to promote the development of uniform filmwise 
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condensation by inducing a stationary vapor condition around the tube. This allowed 
a steam blanket to form around the tube prior to coolant flow initiation. After flow 
initiation the appearance of the condensate film on the tube surface was continuous 
with no breaks. Momentary film instabilities were observed at vacuum conditions at 
pressures below — 20 kPa after a continuous film was established at higher pressures. 
These instabilities seemed to interrupt the film sheet only for an instant and then 
disappear, and may have been caused by vortex shedding of the vapor around the 
tube as already mentioned. A possible mechanism to explain these instabilities is 
that the higher vapor velocity at vacuum conditions momentarily thins the condensate 
film via vortex shedding, yet this thinning effect is overcome by surface tension forces 
in the film sheet which tend to restore the continuous film. These instabilities could 
only be seen for an instant and then would vanish, being very transitory in nature. 
There also did not appear to be any pattern whatsoever to the instability formation. 
The instabilities were not observed at pressures above ~ 20 kPa. 

The data taking regimen for each data set involved starting, then verifying the 
existence of a filmwise condensation condition, then taking data at flow rates (in %) 
of 80, 70, 60, 50, 40, 30, and 20 then back to 80 and 50 to check for repeatability 
within the data set. Two data points were taken at each of the first seven data points 
and one each for the last two, which gave a total of 16 data points. It was usually 
quite clear from the two comparison points, and from data taken previously under 
similar conditions, whether the data set should be rejected or accepted. After tube 


installation, the appearance of one or more small patches (breaks in the film) after 
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— 7-10 hours of operation signaled the beginning of tube contamination which got 


worse with time. The tube would then be removed and cleaned. 


C. TUBES TESTED 

The data taken during this study involved extensive use of the instrumented 
smooth tube (S01) fabricated by Poole [Ref. 6], (six wall thermocouples spaced 60 
degrees apart placed at midwall and midlength). Due to excessive thermocouple 
wear, only 5 thermocouples in this instrumented tube functioned properly. The tube 
was positioned in the apparatus to make optimum use of functioning thermocouples. 
The preferred arrangement placed 4 thermocouples at 10°, 190°, 250°, and 310° from 
the top dead center position of the tube. This arrangement provided readings from 
the top and bottom of the tube and two intermediate points, giving the most accurate 
mean tube wall temperature and best temperature profile readings. At the 
conclusion of this study, the manufacture of the new instrumented tubes was not yet 
complete; progress to date (March 1992) is recorded in Appendix F. 

Data was also taken on a uninstrumented smooth tube (S02) and four finned 
tubes (S03, S04, SOS, S06) with fin spacings of 0.5 mm, 1.0 mm, and 1.5 mm, and 2.0 
mm. All tubes tested were classified as medium tubes, with an outside diameter of 
19.05 mm, and an inside diameter of 12.70 mm. The four finned tubes all had the 
same fin height of 1.0 mm, and fin width of 1.0 mm. Data runs were taken either 
with no insert, with the wire wrap or the Heatex insert installed.By more сїнс 


mixing the coolant an insert significantly increases the inside heat-transfer coefficient. 
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The Heatex insert consists of a central wire core onto which are wound a series 
of wire loops, each inclined at a common angle to the core. The loops come into 
direct contact with the tube wall, and each loop provides a significant amount of 
coolant mixing as the coolant flows through the loop mesh. [Ref. 12]. 

The wire wrap insert was a copper wire spirally wrapped around a central 
stainless steel rod with a uniform pitch. This insert induced a swirling coolant 
motion, which enhanced turbulent mixing within the tube. This particular wire wrap 
insert was the same insert used by Guttendorf [Ref. 3], Coumbes [Ref. 22], and Van 
Petten [Ref. 2] to collect data on the medium family of tubes. 

A summary of data runs is given in Table 2. 


Table 2. DATA RUN SUMMARY 


Tube Filename Pressure Vapor Insert 
(kPa) Velocity 
(m/s) 

501 DSOIMVSH2 11 211 Heatex 
DSOIMASH4 101 1.06 Heatex 
DSOIMASH5 101 1.06 Heatex 
FSOIMASH3 11 2411 Неаїех 

FIMAVSHI 28 3.53 Heatex 
FIMAVSH2 41 2.45 Heatex 
FIMAVSH3 69 1.53 Heatex 
FIMAVSH4 101 1.08 Heatex 
FIMAVSH5 28 1.01 Heatex 
FIMAVSH6 41 1.04 Неаїех 
FIMAVSH7 69 1.02 Heatex 
FIMASW3 101 1.07 wire wrap 
FIMASW4 101 1.07 wire wrap 
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Tube Filename Pressure Vapor Insert 


(kPa) Velocity 
(m/s) 

FIMASWS 101 1.06 wire wrap 
COMPAIWI 101 1.07 wire wrap 
COMPAIW2 101 1.07 wire wrap 

FIMASN4 101 1.08 none 

FIMASNS 101 1.07 none 

FIMASNG 101 1.08 none 
COMPAINI 101 1.08 none 
COMPAIN2 101 1.08 none 

S02 FNMAVSHI 15 6.24 Heatex 
FNMAVSH2 41 2.42 Heatex 
FNMAVSH3 69 LI Heatex 
FNMAVSH4 102 1.03 Heatex 
FNMAVSHS 28 1.02 Heatex 
FNMAVSH6 41 1.03 Heatex 
FNMAVSH7 69 1.02 Heatex 
FNMAVSHS 101 1.07 Heatex 

503 FIMAFOSI 101 1.06 Heatex 

FIMAF052 101 1.06 Heatex 

S04 FIMAF101 101 1.07 Heatex 

FIMAF102 101 1.06 Heatex 

505 FIMAF151 101 1.07 Heatex 

FIMAF152 101 1.07 Heatex 

S06 FIMAF201 101 1.07 Heatex 

FIMAF202 101 1.07 Heatex 
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Tube descriptions: 


S01 
502 
503 
504 
505 


506 


instrumented smooth tube fabricated by Poole 

uninstrumented smooth tube 

0.5 mm finned tube 

1.0 mm finned tube 

1.5 mm finned tube (this tube designated as F096 by Guttendorf, also a 
similar 1.5 mm finned tube designated as F006 by Guttendorf and Van 
Petten was not tested); the 1.5 mm finned tube results were not included 
in subsequent finned tube analysis since the reason for differences in 


experimental results for these two tubes has not yet been resolved. 


2.0 mm finned tube 
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V. THEORETICAL BACKGROUND AND DATA REDUCTION PROCEDURES 


A. THEORETICAL BACKGROUND 


The overall or total resistance to heat transfer from vapor to coolant consists 
of the sum of the vapor side resistance (R,), the tube wall resistance (R,), and the 


coolant side resistance (R,); this neglects any fouling resistance since clean tubes are 


always used. 
Ra =, +R, R, (5.1) 


The vapor and coolant side resistances are convective in nature and may be 
expressed by the reciprocal of their respective heat-transfer coefficient and surface 


area product. 


1 
К = —— (5.2) 
° h,A, 
1 
I aS (5.3) 
һ,А, 
where: 
R, = outside vapor side resistance to heat transfer (K/W) 
h, = outside heat-transfer coefficient (W/m^K) 
A, = effective outside surface area (m^) 
R; = inside coolant side resistance to heat transfer (K/W) 
h; = inside heat-transfer coefficient (W/m^K) 
A;  - effective inside surface area (m^) 


The tube wall resistance is conductive in nature and is represented by the radial 


conduction equation. 
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р 
In— 
D 





| (5.4) 
R= 
271, 
where: 
К, = tube wall resistance (K/W) 
D, = outside or root (if finned) diameter (m) 
D; = inside tube diameter (m) 
L = active condensing length (133 mm 
К, - thermal conductivity of tube wall (W/m-K) 
The effective outside area of the tube is calculated using the following 
expression: 
A,="D_L (5.5) 
where: 
D, = outside or root (if finned) diameter (m) 
L = active condensing length (133 mm) 


The effective inside area includes the inside surface area involving the active 
condensing length and the inside surface area of the insulated inlet and outlet 
portions of the tube. These portions of the tube act as fins and remove heat via 
axial conduction. The extended fin assumption with associated fin efficiencies was 


used to account for these end losses. 


A; 2 1 DL * n, L, *n, L,) (5.6) 
where: 
D,  - inside diameter of tube (m) 
Г = active tube condensing length > 
7, = fin efficiency of inlet portion of tube 
Г. > length of inlet portion of tube (m) 
7, == fin efficiency of outlet portion of tube 


length of outlet portion of tube (m) 


The overall thermal resistance to heat transfer may be expressed by: 
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=R,+R,+R, (5.7) 


о о 


Substituting equations (5.2) and (5.3) into equation (5.7) yields: 











! = 1 * К, T 1 (5.8) 
U,A, h,A, һ,А, 
where: 
U, = overall heat-transfer coefficient (W/m К) 
A, = effective outside surface area ( m) 
h, = outside heat-transfer Coefficient (W/m^K) 
R, = Tube wall resistance (K/W 
h; = inside heat-transfer coefficient (W/m^K) 


The single tube condenser apparatus uses the log mean temperature difference 


(LMTD) analysis for calculation of the heat transfer between the hot vapor and cold 


coolant. 
Q=U,A,(LMTD) (5.9) 
where: 
Q = heat transfer rate to the cooling water DR)? 
U, = overall heat-transfer coefficient (W/m* 
A, = effective outside surface area (mî) 


LMTD = log mean temperature Ta between vapor and coolant (K) 


The log mean temperature difference (LMTD) is given by: 
(7; - T) 
In I. 2 Т, 
Tsa 1 


LMTD = Ж” 
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where: 


T, = coolant inlet temperature (K 
T, = coolant outlet temperature (K) 
Ты = vapor saturation temperature (K) 


In this and previous studies at NPS, the quartz thermometer output for T, and 
T, were used in the calculations for the coolant temperature rise, and the saturation 
temperature, T,» was measured using the vapor thermocouple (channel 40). 

The total heat transfer across the tube is experimentally determined by 
measuring the mass flow rate of fluid through the tube and its accompanying 


temperature rise. 


Q-mc,(T, - T) (5.11) 
where: 
Q = heat transfer rate (W) 
m = mass flow rate of coolant (kg/s) 
C, = specific heat of coolant at constant pressure (J/kg: K) 
T, = coolant inlet temperature ( E 
T, = coolant outlet temperature (K) 


Equation (5.11) may be used directly from the experimental data. The 
resultant heat transfer rate, Q, is then substituted into equation (5.9) to find the 


overall heat-transfer coefficient, U,. 


D — (5.12) 
A,(LMTD) 
where: 
Q = heat transfer rate from eq. (5. 11) (W) 
A, = effective outside surface area (m^) 
LMTD = log mean temperature difference; eq. (5.10) (K) 
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Since R,, U,, A,, and A; are known quantities, this leaves only two unknowns 
in equation (5.8), h, and h,, the outside and inside heat transfer coefficients. 

Often the coolant side thermal resistance is dominant, and inserts such as those 
mentioned in section IV C are used to lower the inside resistance. This allows a 
more accurate computation of the outside heat-transfer coefficient, h,, when using 
the modified Wilson plot technique mentioned in section V B. Vapor side heat 
transfer may also be enhanced through the use of fins, drainage strips, or dropwise 


condensation promoters. 


B. MODIFIED WILSON PLOT TECHNIQUE 

The ideal way to solve for h, and h; in equation (5.8) is through the use of 
instrumented, which accurately determine a mean tube wall temperature. The inside 
and outside mean tube wall temperatures may then be obtained directly by assuming 
a linear temperature profile across the wall. Since the vapor temperature and mean 
coolant temperature are known, the inside and outside heat-transfer coefficients may 


then be calculated directly using equation (5.13). 


q-hAT (5.13) 
where: 
q = heat flux (W/m?) 
h = heat-transfer coefficient (h; for inside, h, for outside) (W/m^K) 
AT = temperature difference across resistive medium (AT = T,,-Tyanousside 


for the vapor side, and AT = T, Toons for the coolant side) 


allinside ^ *co 


For data collection on a large number of tubes, the use of instrumented tubes 
Is impractical due to the high cost and difficulty involved in manufacturing so. many 
tubes. Therefore, the modified Wilson plot technique was developed, which solves 


for the inside and outside heat-transfer coefficients simultaneously without using wall 
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thermocouples. To obtain the most accurate results with this method, it is necessary 
that the inside and outside coefficients be relatively equal in magnitude. 

The modified Wilson plot technique requires that the "form" of the equation 
for both the inside and outside heat-transfer coefficients be known. The Nusselt 
theory and Sieder-Tate correlation are used to represent the "form" of the outside 
and inside heat-transfer coefficients respectively. The Nusselt theory when based on 


q can be represented by: 


Ер аһ 
Л, =a 22278 =aZ (5.14) 
p,D,q 

or 
where: 

h, = outside heat-transfer coefficient (W/m*K); based on q 

a = dimensionless coefficient 

К, = thermal conductivity of condensate film (W/m-K) 

р; = condensate film density (kg/m ), 

Р = gravitational constant (9.81 m/s’) 

g = Specific enthalpy of vaporization (J/kg) 

и; = dynamic viscosity of condensate film (N-s/m’) 

D, — outside tube diameter (m) 

q = heat flux based on outside area (Q/A,) (W/m’) 

The Sieder-Tate correlation may be represented by: 

k u yn 
h, - C, — Re°* Pr T - CQ (5.15) 
D i Hy 

where 

h; = inside heat transfer coeffieient (W/m^K) 

C; = Sieder-Tate leading coefficient 

К. = thermal conductivity of cooling water (W/m:-K) 

D, = inside tube diameter (m) 

Re = Reynolds number 

Pr = Prandtl number 

и. = dynamic viscosity of cooling water at bulk temperature (N-s/m’) 


4] 


dynamic viscosity of cooling water at mean inner wall temperature 
(N:s/m^) 
substituting equations (5.14) and (5.15) into equation (5.8) gives: 


Hw BE 




















AZ 
— үен (5.16) 
© ОА, « 
By letting: 
y-|4 -R AZ (5.17) 
X= A (5.18) 
А0 
— (5.19) 
С, 
1 
b=— (5.20) 
а 
then equation (5.16) reduces to: 
Y=mX+b (5.21) 


The parameters () and Z are both temperature dependent, and must be 


determined iteratively. A least-squares fit of equation (5.21) is used to determine 


C; and a. Once C; is known, h, can be calculated using equation (5.15). With h; and 
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U, known, the value of h, can then be easily determined by rearranging equation 


(5.7) as in equation (5.22), or by using equation (5.14) with the value of a known. 


(5.22) 








C. INSTRUMENTED TUBE IMPROVEMENTS FOR DATA REDUCTION 

In previous work at NPS, the standard form of the Sieder-Tate equation was 
used with a Reynolds number exponent of 0.8, as in equation (2.10). One of the 
aims of this thesis was to use an instrumented tube to directly determine the inside 


and outside coefficients, h; and h, , and then use the data to determine a more 


о › 
"exact" form of the Sieder-Tate-type equation to be used for each insert. The 
coolant side correlations mentioned in section II C were based on a long, straight 
entrance length. The sharp 90° bend just prior to the test section tube entrance 
undoubtedly creates entrance effects which lead to discrepancies between our 
experimental data and heat transfer behavior predicted by the Sieder-Tate 
correlation. 


Assuming the final form of the inside heat-transfer correlation to be that of 


equation (2.10) gives the following: 


7044 
Nu - C, Re" Pr | (5.23) 
. Hy 
where: 
m = Reynolds number exponent to be determined 
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Rearranging equation (5.23) gives: 


Nu = С Ке" 


p, in( Hc 014 o! (5.24) 
Bu 


Taking the natural log of equation (5.24) yields: 


In | ——————_—_ | = In Ci+ m In Re 
Ч u Г (5.25) 
Prizes 


Equation (5.25) is in the form of a linear equation, and by plotting /n 
(Nu/(Pr'^(uJu,)"*)) versus In Re, the slope and intercept, namely the Reynolds 
number exponent and Sieder-Tate coefficient may be determined from the 
instrumented data. 

With the unknown parameters of equation (5.23) determined, the new inside 
heat-transfer correlations (one for each insert) could be used in the data reduction 
program to give the value of the inside heat-transfer coefficient, h,, directly and to 


provide a more accurate calculation of the outside heat-transfer coefficient h,. 


D. ENHANCEMENT RATIO 
Following the development of Van Petten [Ref.2], and Nusselt theory, 
experimental data can be curve fitted, using a least-squares analysis, to an equation 


of the following form: 
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а = аАТ” (5.26) 


where: 


AT 


Substituting equation (5.12) into his expression yields: 


heat flux based on outside area (Q/A,) гт 
temperature drop across the condensate film (K) 


Л, -аАТ"! (5.27) 


From Nusselt theory п - 0.75, therefore the enhancement ratio, based on 


constant temperature drop across the condensate film, can be expressed as: 


h а 
Е ој 
Єт = — (5.28) 
05 4, 
where: 
€r = enhancement ratio based on constant temperature drop across the 
condensate film 
f = subscript denoting finned tube 
5 — subscript denoting smooth tube 
h,  - outside heat-transfer coefficient (W/m^K 
a — constant of proportionality introduced in equation (5.26) 


Also for constant AT; using equation (5.14): 

1 1 

“ДААЛТ ала, 
а|адТ &,| а; 


If the heat flux is kept constant the values of Z, and Z, remain equal, which 


q. 1/3 


(5.29) 











results in equation (5.30). 
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Combining equations (5.29 and (5.30) gives the relationship between €,, and 
e, in equation (5.33). 











1/3 

“шинэ МЭ (5.31) 
En 
с. а, а, 
[3 

c e Ref (5.32) 

4 o a, 
€= GO (5.33) 


Note that €,; and е, аге independent of q and AT. 
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VI. RESULTS AND DISCUSSION 


A. DROPWISE CONDENSATION 

As mentioned in section IV B, a dropwise condensation condition was obtained 
initially. Data were taken during dropwise conditions at vacuum and atmospheric 
pressure with the instrumented tube, and this data is compared to filmwise data in 
Figure 4. The figure shows a marked contrast between filmwise and dropwise 
condensation data. The dropwise heat-transfer enhancement, compared to the 
filmwise data, varied from ~2 to 7 for vacuum conditions and ~9 to 10 for 
atmospheric conditions. Marto et al [Ref. 23] studied the use of organic coatings for 
the promotion of dropwise condensation of steam, and obtained an outside heat- 
transfer coefficient of —55 kW/m^K for a Fluoroacrylic coating compared to —30 to 
85 kW/m^K found under vacuum conditions in this study (both taken at P— 11 kPa 
and vapor velocity —2 m/s). The difficulty in accurately measuring the temperature 
drop across the condensate film for the dropwise condition is illustrated by the large 
amount of scatter in the dropwise data; this problem is caused by the lack of a stable 
film and the intermittent presence of drops near the instrumented tube wall 
thermocouples. The exact cause of the dropwise conditions was never determined. 
Clearly, some organic contamination either from the boiler feed water or from the 
gasket material was depositing on the test tube. Since this thesis was devoted to 
filmwise condensation, great efforts were made to clean the test tube and prepare 


it chemically so that the condensate would wet the surface. 
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Figure 4. Comparison of Dropwise and Filmwise Condensation Data 
Instrumented Tube, Heatex Insert) 
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B. INSTRUMENTED TUBE RESULTS 

After resorting to the oxide coating procedure mentioned previously to mitigate 
the dropwise contamination problem, a series of runs were made with the 
instrumented tube fabricated by Poole [Ref. 6] with the Heatex insert installed. 
Using a mean wall temperature, the inside and outside heat-transfer coefficients 
could be evaluated directly. 

Figure 5 shows filmwise condensation data taken at various pressures and vapor 
velocities. The four data sets for increasing vapor velocity (from 1 to 3.5 m/s) and 
decreasing vapor pressure (from 101 to 28 kPa) show the effect of vapor shear 
thinning the condensate film, giving an increase in the outside heat-transfer 
coefficient. This series of data runs was taken by maintaining the heater voltage at 
175 volts and adjusting coolant flow through the auxiliary condenser to control 
pressure. 

The four data sets in Figure 5 with constant vapor velocity (—1 m/s) and 
decreasing vapor pressure (from 101 to 28 kPa) show reduction of the outside heat- 
transfer coefficient with decreasing saturation temperature (due to decreasing 
saturation pressure). This effect is thought to be due to an increase in condensate 
viscosity at lower temperatures which tends to prevent the condensate from flowing 
around and draining from the tube as easily as at higher temperatures. The resulting 
condensate film thickening provides an additional resistance to heat transfer thereby 
lowering the outside heat-transfer coefficient. This series of data runs was taken by 
both adjusting heater power and auxiliary condenser coolant flow to obtain the same 


pressure conditions as above with a constant vapor velocity of —1 m/s. 
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50 





Figure 6 shows how increased pressure, with vapor velocity held constant, gives 
a larger increase in the outside heat-transfer coefficient over the predicted Nusselt 
value. This discrepancy is also thought to be caused by the viscosity effect since the 
Nusselt treatment seemingly takes into account the other possible causes. Also 
depicted more clearly is the dramatic enhancement due to vapor shear effects for 
increased vapor velocity at constant pressure. 

As mentioned previously, the use of inserts allows more effective mixing of 
coolant and facilitates greater accuracy in the calculation of the outside heat-transfer 
coefficient. Figure 7 shows the outside heat-transfer coefficient determined from the 
use of the two inserts together with no insert. The results for the wire wrap and 
Heatex inserts are closely grouped, whereas the no insert data shows much greater 
scatter. This seems to indicate that the use of an effective insert does indeed 
enhance accuracy for the instrumented tube data. 

The mean temperature difference across the condensate film for the no insert 
case has a significantly lower value than either of the insert cases; with the outside 
and tube wall conditions unchanged, this indicates a higher inside resistance for the 
no insert case. To illustrate this point a mid-range instrumented tube data point at 
the same conditions for each insert case is shown in Table 3 (P,,—101 kPa, 


V...,-~ 1.1 m/s, V 


I ~2.75 m/s). It shows that both inserts have a comparable effect 


coolant 
(slightly better with Heatex) and roughly provide a factor of two enhancement over 
the no insert case in the inside heat transfer coefficient. Due to the increased inside 


resistance for the no insert case, the heat flux shows a decrease of about 1546. 
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Figure 7. Comparison of Steam Heat-Transfer Coefficients at Atmospheric 
Conditions (P,,,— 101 kPa, У, — 1.1 m/s) for Three Insert Conditions 
(Smooth Instrumented Tube) 


53 


Table 3. PERFORMANCE COMPARISON OF TUBE INSERTS 
AT ATMOSPHERIC CONDITIONS 


Heatex Insert Wire Wrap Insert No Insert 
h; (KW/m^K) 30.85 27.63 15.31 


q (kW/m?) 481.4 477.7 409.6 
h, (kW/m?-K) 9.86 9.82 10.75 
AT(K 48.82 48.65 38.10 





Figure 8 shows the wall temperature profiles for the data points listed in Table 


3 along with a vacuum data set (P,,—28 kPa, V,,,,,—1 m/s, V ~2.75 m/s) for 


vapor coolant 


comparison, where 0 degrees is at the top dead center position of the tube. The 
shape of the temperature profiles shows the effect of condensate film thickening 
toward the bottom of the tube. The higher resistance through a thicker condensate 
film results in a lower tube wall temperature toward the bottom of the tube as 
shown. 

As expected, the wire wrap insert and Heatex insert temperature profiles at 
atmospheric pressure are very similar. The no insert profile has a shape similar to 
the wire wrap and Heatex profiles, but the mean wall temperature is about 11 K 
higher and shows the effect of the increased inside resistance to heat transfer. The 
vacuum run temperature profile shows the effect of a lower temperature gradient 
between the steam and coolant; the lower heat-transfer potential results in lower 
heat fluxes and a flatter temperature profile. 

Figure 9 shows the comparison of the instrumented data with the predictions 
of Nusselt, Fujii, and Shekriladze and Gomelauri covered in section II B. The data 
depicted ranges from P,,, — 28 kPa and V,,,,, — 3.5 m/s to P,,, — 101 kPa and V, 


vapor apor 


= ].1 m/s. The data seems to follow the Shekriladze and Gomelauri prediction the 


closest, but is also very near the Fujii prediction. 
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Figure 9. Comparison of Experimental Results with the Predictions 
of Nusselt, Fujii, and Shekriladze-Gomelauri. 
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C. INSIDE HEAT-TRANSFER CORRELATIONS FROM INSTRUMENTED 

TUBE RESULTS 

Тһе instrumented tube provides a method of determining an inside heat- 
transfer coefficient for each insert from direct wall temperature measurements. 
Wanniarachchi et al [Ref. 8], in a similar effort to resolve the differences between 
the Sieder-Tate correlation and his experimental results, developed a correlation 
based on the Sieder-Tate expression. Using an intercept form, with data taken on 
the same apparatus as used in this study, and a least-squares fitted leading 


coefficient, the correlation took the following form: 


0.14 
Nu =0.064 Re Pr m +26.4 (6.1) 
Bw 


However, Rouk [Ref. 4], using an optimization technique, showed that the value of 
the intercept had little effect on the results, and it was the accurate determination 
of the Reynolds exponent that was more critical. 

With a view to finding the appropriate Reynolds exponent and leading 
coefficient for equation (5.23), the data for no insert, wire wrap, and Heatex insert 
were plotted as /n Re versus /n (Nu/Pr'" (i1, )^ ^) as explained in section V C. The 
plotted data are shown in Figures 10, 11 and 12 for no insert, wire wrap insert, and 
Heatex insert respectively. A line of best fit (typically with a regression coefficient 
of 0.99) was used to obtain the value of the intercept. Figure 13 shows the plot of 
all three cases on the same graph for comparison. Note that the insert data is 
closely grouped on the upper regions of the graph when compared to the no insert 
data. The increase in the Nusselt number again indicates more efficient inside heat 


transfer for the insert vice the no insert case. The Reynolds exponent, or slope, 
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Figure 10. Log-Log Plot of Re versus Nu/Pr'^(1,/u, ^ for No 
Insert (Smooth Instrumented Tube 
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differs from the Sieder-Tate-type equation of 0.8 in all three cases. The following 
derived correlations apply specifically only to the medium tube, but should be 
applicable to any tube with the same inside diameter. 


The no insert inside heat-transfer correlation had the form: 


0.14 
Nu 20.013 Re? 9 pr "s (6.2) 
By 


Equation (6.2) was used to reprocess some current smooth tube runs, with no insert, 
as well as those of Guttendorf [Ref. 3] and Van Petten [Ref. 2] to check their values 
for the smooth tube Nusselt coefficient a. 


The wire wrap insert inside heat-transfer correlation had the form: 


0.14 
Nu -0.052 Re99€ p; e (6.3) 
Hy 


Equation (6.3) was used to reprocess previous data taken on the same apparatus by 
Guttendorf and Van Petten, who used the wire wrap insert for their finned and 
smooth tube experiments, with a view to checking their results with this new 
correlation. Table 4 shows a comparison between the new wire wrap insert 
correlation (eq. 6.3), and Wanniarachchi's correlation (eq. 6.1) both with and without 
the intercept value included. This comparison shows that equation 6.1 is actually 
more accurate without the intercept value included when compared to the results of 


the new wire wrap insert correlation. 
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Table 4. COMPARISON OF EQUATIONS (6.1) AND (6.3) 
FOR Pr'” ^2 1.4 (held constant for comparison 


Reynolds Number 20,000 | 30,000 | 40,000 | 50,000 

Nusselt Number 2757 368.4 456.9 541.0 
Ми=0.064 Re?*Pr'^ J^ - 26.4 

Nusselt Number 247.3 342.0 430.5 514.6 
Nu=0.064 Re®®Pr'” Л” 
244.9 341.5 432.3 519.1 


The Heatex insert inside heat-transfer correlation had the form: 
















Nusselt Number 
Nu=0.052 Re?*?pr!? 


0.14 
Ми = 0.22 пери < (6.4) 
В» 


Equation (6.4) was used to reprocess all Heatex data. 

Memory [Ref. 12] conducted condensation experiments on a different apparatus 
with a smooth instrumented tube. He also determined the Reynolds number 
exponent for no insert, wire wrap insert (made locally and somewhat different from 
the one used in this study), and Heatex insert. The exponents he obtained are 
reported in Table 5. The Heatex insert exponent of 0.68 compares very favorably 
with the value of 0.69 found in this study. The wire wrap exponent of 0.73 was well 
below the value of 0.82, most likely due to differences in the insert. The no insert 
case gave a value of 0.85 and compared well to the value of 0.89 found in this study. 
The difference in the no insert exponent is most likely due to the difference in the 
tube entrance region; Memory used a long run of straight pipe for the tube entrance, 
whereas the present work had a sharp bend just prior to the condenser tube. 

Rouk [Ref. 4] used the instrumented smooth tube data of Georgiadis [Ref. 5] 


to the find the appropriate value of the Reynolds number exponent for the no insert 
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and wire wrap cases. His results compare quite well with this study and are also 


given in Table 5. 


ANL (Argonne National Laboratory) [Ref. 21] conducted an assessment of 


heat-transfer correlations for turbulent pipe flow with water to determine the best 


correlation(s) on which to base their design of Ocean Thermal Energy Conversion 


(OTEC) heat exchangers. ANL used two shell-and-tube heat exchangers, with no 


inserts, for analysis and reported the following: 


where: 


where: 


AA 


кә 


The Dittus-Boelter (eg 2.7) and Sieder-Tate (e 
predicted the data by 
design. 


1 correlations under- 
ere 


o to 1596 and were consi too conservative for 


Overall, the "best" correlations were found to be Petukhov-Popov (eq 6.5) and 
Sleicher-Rouse (eq 6.6), both of which showed excellent agreement (+ 5%) 
with the experimental data (at Pr=6.0 and Pr=11.6). 


Nus. CIO ROB ___ (6.5) 
K, * K,(e[8) ^ (Pr?? -1) 
(valid for 0.5 < Pr < 2000 and 10° < Re < 5x106) 
я log, Re - 1.64)? 
+ 3.4є 
11.7 + 1.8 PES 
Nu =5 +0.015 Ref Pr, (6.6) 


(valid for 0.1 < Pr < 10° and 10° < Re < 10°) 


0.88 - 0.24/(4+ Pr,) 
1/3 OSE 
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3. The most accurate correlations (i.e. Petukhov-Popov and Sleicher-Rouse) 
seem to yield effective Reynolds exponents in the neighborhood of 0.85 
(uncertainty range: m = 0.82 to 0.88). 


4. The potential sources of uncertainty in the Wilson procedure included 
waterside flow maldistribution, entrance effects, experimental error in U,, and 
the uncertainty in the Reynolds number exponent. Of these, they concluded 
that the uncertainty in the Reynolds number exponent was, by far, the most 
significant. In fact, the results of the Wilson procedure were found to be 
highly sensitive to the value of the Reynolds number exponent. 


Table 5. COMPARISON OF REYNOLDS NUMBER EXPONENTS FOR 
SIEDER-TATE-TYPE CORRELATIONS 


Experimental Rouk 
Data 


No Insert 0.89 0.90 


Wire Wrap 0.82 0.78 
Heatex 0.69 --- 





D. ANALYSIS OF SMOOTH TUBE RESULTS 

When using the modified Wilson plot technique to reprocess data files, the 
solution option can be specified to use either the stored value of the Sieder-Tate 
coefficient (for direct computation of h,) or let the coefficient value "float", which 
allows the program to calculate its own value of the coefficient. In order to 
determine which method was most accurate, the instrumented data files were 
reprocessed using each method and then compared with the values of the heat- 
transfer coefficient which were obtained by direct measurement of the tube wall 
temperature. A high, medium, and low coolant flow rate was chosen from each run 
to facilitate the comparison. The results were tabulated and are shown in Appendix 
D; it can be seen that the fixed coefficient method yielded the more accurate results 
at least 75% of the time. The mean error of the fixed method was «2.0%, and that 


of the floating method was +5.4%. The error for the lowest coolant flow rate (Re 
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<20,000) was noticeably higher than for higher coolant flow rates for both methods. 
The choice of using the fixed coefficient method represents a departure from the 
practice of previous researchers on this apparatus who exclusively used the floating 
coefficient method. 

Prior to the instrumented tube runs reported in Figure 5 (from which the new 
correlations were empirically derived) a series of data runs were made using a plain 
smooth uninstrumented tube (S02) of the same dimensions using the Heatex insert. 
The plain smooth tube data was then reprocessed using the new correlation for the 
Heatex insert and the results are plotted in Figure 14. These data sets were taken 
at the same conditions as the instrumented data of Figure 5 except for the set at the 
highest vapor velocity of 6.2 m/s vice 3.5 m/s for the instrumented tube. Similar 
effects of vapor shear and vapor pressure, as mentioned previously for Figure 5, are 
clearly seen, and again illustrate the vapor shear effect on the outside heat-transfer 
coefficient. 

With the exception of the two data runs at high vapor velocity, the data from 
Figures 5 and 14 are shown together in Figure 15. The close agreement of the 
reprocessed plain smooth tube data with the instrumented smooth tube data allows 
a high degree of confidence in the accuracy of the new correlations and the choice 
of the fixed coefficient method. 

To provide a baseline from which to evaluate finned tube performance it was 
necessary to obtain the smooth tube Nusselt coefficient, a, for the specific conditions 
under which the comparison was to be made. The condition chosen was atmospheric 
pressure (101 kPa) and a vapor velocity of ~1 m/s. 

For 8 complete sets of data the average value of a, using the fixed ШЕШІ was 


found to be 0.876. The average value for each data set was found by taking the 
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reprocessed plain smooth tube data 
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Figure 15. Comparison of Instrumented Smooth Tube Results with Non- 
Instrumented Smooth Tube Data After Reprocessing with the 
New Heatex Insert Inside Heat-Transfer Correlation. 
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(K) 


( Tsat - Twall ) 


measured value of h, for each data point, dividing it by the Nusselt theory prediction 
of h,, and then multiplying by 0.728 (the Nusselt coefficient); the average of all the 
data points in the set was then taken. Interestingly the Wilson plot floating 
coefficient method gave a value of 0.835, somewhat lower than the average value. 
Originally it was thought that this discrepancy might be due to "outlier" data points 
(high or low coolant velocity) in each set. However, removing the highest or lowest 
coolant flow rates within a set had little or no effect on the Wilson plot result. The 
reason for the discrepancy is still not known and merits future study. 

The value of «a was calculated for several other flow conditions; these are 
shown in Table 6. The trend of the readings, like that of Figure 6, shows that vapor 
velocity has a much greater effect than pressure on the value of a, as expected. 


Table 6. SMOOTH TUBE a SUMMARY; EFFECT OF PRESSURE AND 
VAPOR VELOCITY 


File Name Vapor Velocity 
(m/s) 


FIMAVSH1 3.5 
FIMAVSH2 25 
FIMAVSH3 155 


FIMAVSH4 Ыы 
FIMAVSH7 1.0 
FIMAVSH6 1.0 
FIMAVSHS 1.0 





E. ANALYSIS OF FINNED TUBE RESULTS 

With an accurate value of a, and the newly determined inside heat-transfer 
correlations, the medium family finned tube data of Van Petten [Ref. 2] was 
evaluated. Figure 16 shows the data Van Petten reported in his thesis; it also shows 


his data after being reprocessed using the new wire wrap insert correlation with 
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Figure 16. Comparison of the Steam Heat-Transfer Enhancement Data of Van 
Petten, for the Medium Finned Tube Family, Using the Modified 
Wilson Plot and New Wire Wrap Insert Inside Correlation. 
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both the fixed coefficient and Wilson plot floating coefficient methods. Since Van 
Petten used the Wilson plot method, it is not surprising that the original thesis data 
and new Wilson plot floating coefficient data are comparable since the Reynolds 
exponent only varied from 0.8 to 0.82. The fixed coefficient method enhancement 
is substantially higher than the Wilson plot results. 

Since the assertion is that the fixed coefficient method is more accurate than 
the Wilson plot method, then the conclusion must be that the enhancement for this 
set of finned tubes is actually higher than previously reported. 

During this study, limited medium finned tube experiments were conducted for 
purposes of comparison. Figure 17 shows the comparison between this data and the 
newly reprocessed data of Van Petten (using the fixed coefficient method) and shows 
reasonable agreement. To more clearly illustrate this point, the data taken on the 
2.0 mm fin spacing tube has been given in more detail in Figure 18. Excellent 
agreement is seen between the experimental results of this study and that of Van 
Petten using the known inside heat-transfer correlation with the fixed coefficient 
method. Again, as shown in Figure 16, the Wilson plot prediction is significantly 


below the fixed coefficient results. 
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Figure 17. Comparison of the Steam Heat-Transfer Enhancement Data 
of Van Petten and Swensen for the Medium Finned Tube Family. 
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Figure 18. Comparison of the Steam Heat Transfer Data of Van Petten and 
Swensen for the 2.0 mm Fin Spacing Medium Tube. 


73 


K ) 


( 


( Tsat - Twall ) 


А. 


VII. CONCLUSIONS AND RECOMMENDATIONS 


CONCLUSIONS 


The inside heat-transfer correlation is highly sensitive to the Reynolds number 
exponent. 


Each insert condition must be analyzed separately to determine the 
appropriate "form" of the inside heat-transfer correlation. 


Calculations based on a known inside heat-transfer correlation are more 
accurate than modified Wilson plot results. 


Armed with accurate inside heat-transfer correlations, previous data may be 
reprocessed to give more accurate results. 


The source of contamination in the test apparatus, which has caused a 
dropwise condensation problem for a number of years, is most probably due 
to a contaminated distilled water source. 


RECOMMENDATIONS 


Reprocess all previous medium and large diameter finned tube data using the 
fixed coefficient method to obtain more accurate results. 


Continue with construction of smooth instrumented tubes of different 
diameters (1.e. small, medium, and large) to confirm the medium tube results 
and develop correlations specifically for the small and large diameter tubes. 


Construct one representative instrumented finned tube to test the validity of 
applying instrumented smooth tube results to finned tube data. 


Test representative water samples that have been collected from both the 
distiller and boiler to confirm the presence of impurities and validate their 
origin. 


Replace current distiller with a deionized pure water source (either 
commercial purchase of water or new distilling apparatus). 
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APPENDIX A. PHYSICAL AND THERMODYNAMIC PROPERTIES OF 
WATER 


The physical and thermodynamic properties of water were based on the 


following equations: 
Heat Capacity (J/kg K): 
C,-4211 - T «[2.268x 10?- T « (4.424x10? -2.714x10 7 « T)] (A.1) 


where: T = temperature (celsius) 


Dynamic viscosity, (kg/m:s): 


1 =(2.4x 10°5) + 101247.37+ 133.19] (A.2) 


Thermal conductivity (W/m:K): 
k = -0.9225 + x «(2.8395 - x «(1.8007 -x «(0.5258 - 0.0734 «х))) (A.3) 


where: x = (T+273.15)/273.15 


Density (kg/m): 


p 2 999.5295 « T« (0.0127 -T* (5.4825x10 ? -T« 1.2341x10 5)) (A.4) 


> 


Latent heat of vaporization (J/kg): 


h,, = 2477200 - 2450 « (T-10) (А.5) 


Fluid Enthalpy (J/kg): 


h,= T* (4.2038 - T» (5.8813x10 *- T*4.5516x10 9) (A.6) 


saturation pressure (Pa): 


P. = 22120000 x Pr (A.7) 


where: 
Pr = е” 


Br = SUM /[T, (1+ 4.1671 (1-T,) + 20.9751 (1-T,)’) - (1-T,)/ 
(01.0 х 10%(1-Т,) + 6)] 


50М - сон 26.0802 (1-Т): - 168.1707 (1-T,)’ + 64.2329 (1-T,)* - 
6 (1- r) 


118.964 


T, = (T+273.15) / 647.3 


Ї 
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APPENDIX B. SYSTEM CALIBRATIONS AND CORRECTIONS 
B.1 Thermocouple and Quartz Thermometer Calibration 
Several different thermocouples and the quartz crystal thermometer (HP 
2804A, Ser. No. 2244AD1192) were calibrated against a platinum resistance probe 
in a mixed isothermal ethylene glycol bath with a Rosemont Galvanometer model 


920A commutating bridge (Ser. no. 013494) from 23 to 26 September, 199]. 


NOMENCLATURE 

T1 Quartz Thermometer Measuring Probe T1 (Ser. No. 2120A-00707, with dial 
setting 481) 

T2 Quartz Thermometer Measuring Probe T2 (Ser. No. 21204-60459 with dial 
setting 510) 


Т1-Т2 The T1-T2 reading on the quartz thermometer 

1-55 Large diameter metal sheath thermocouple (diameter — 0.040") 

T-56 Small diameter metal sheath thermocouple (diameter = 0.020") 

T-57 Old Thermocouple (taken from rig during disassembly set to HP 3497A 
internal zero); (Type: Omega, TT-T-30, Lot# HCP093HC0306) 

T-58 New I Thermocouple (set to HP 3497A internal zero); (Type: Omega, TT- 
T-30 SLE, Lot OCP1453PTCC01473P) 


T-59 New 2 Thermocouple (referenced to ice bath zero) 

T-60 10 Junction Thermopile #1 (referenced to ice bath zero) 
T-61 10 Junction Thermopile #2 (referenced to ice bath zero) 
T-62 10 Junction Thermopile #3 (referenced to ice bath zero) 


The calibration was performed by taking the bath temperature up from 290 K 


to 393 K by 5 K increments then back down to check for hysteresis; no hysteresis was 


TH 


observed. The data was fitted to a fifth-order polynomial (regression coefficient — 


1.000 for each polynomial fit) in each case with the following results: 


T55 = 273.15 +(2.5943e -2)V-(7.2671-7)V? 
+ (3.2941 -11)V3-(9.7119-16)V* (В.1) 
+ (9.7121 -20)V? 


156 - 273.15«(2.5878e -2)V-(5.9853e -7)V? 
- (3.1242e -11)) V? «(1.3275e - 14) V^ (B.2) 
- (1.0188e-18)V? 


T57 = 273.15 +(2.5923e-2)V-(7.3933e -7)V? 
+ (2.8625e -11)V° +(1.9717e-15)V* 
- (2.2486e-19)V? 


(B.3) 


158 - 273.15«(2.5931e-2)V-(7.5323e -7D)V? 
+ (4.0567e -11)V?-(1.2791e -15)V^ 
* (6.4402e -20)V? 


(B.4) 


T59 - 273.15*(2.5471e -2)V-(3.7621e -7)V? 
- (1.0105e -10)V? «(2.3928e -14) V^ 
- (1.6440e - 10) V? 


(B.5) 


760 = 273.17+(2.5571e-2)V-(1.9980e -7)V? 
- (1.6385e -10)V? «(2.6164e -14)V* (B.6) 
- (1.0295e-18)V? 


T61 » 273.15*(2.6119e -2)V -(9.0449e - 7) V? 
* (1.1214e- 10)V? -(1.5623e -14)V* 
* (1.0646e -18)V? 


(B.7) 
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T62 - 273.15*(2.5996e -2)V -(7.4405e - T) V? 
* (2.4733e-11D)V? «(3.32366 - 15)V* (B.8) 
- (3.7460e -19)? 


For 755, T56, T57, T59; V— Voltage in microvolts (eg. 0.010 volts = 10,000 
microvolts). 

For T60, T61, T62; V= Voltage in microvolts/10 (eg. 100,000 microvolts + 10 
= 10,000). 

The HP 2804A Quartz Crystal Thermometer was also calibrated on 26 


September, 1991 with the results summarized in Table B.1. 


Table B.1. QUARTZ THERMOMETER CALIBRATION DATA 


Reference Ti T2 Error 
Temperature (deg C) (dec C) 
(deg C) 


Ш979 18.003 0.022 
18.977 18.997 0.021 
19.984 20.001 0.018 


20.980 21.007 0.028 
22.936 22 959 0.025 
25.014 25.032 0.023 
30.140 30.153 0.024 
34.977 34.983 0.021 





The mean error for the Quartz Thermometer as shown in the Table B.1 is 
~ +0.02 deg C. However, it is well noted that T7-72, the critical measurement, 
indicates an apparent error of less than +0.005 deg C when the temperature reading 
is near or below 25 deg C; T] and 72 tracking well together lowers the error 


estimate. 
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B.2 Flow Meter Calibration 

The flow meter calibration for coolant flow through the single horizontal tube 
was completed on Oct 28, 1991 using a stop watch, portable tank, and a Toledo 
model 31-0851 IV, Se. No. 1326 scale with 1/10 pound graduations. The following 


relation was obtained via linear regression: 


m = 6.7409F + 13.027 
B.9 
where: 


= mass flow rate oo per second) 
= flow meter reading (eg. 10%=>10) 


m 
F 
The applicable range of the calibration was 10% to 95%. The water 
temperature on the day of the calibration was 17.5°C (290.6 K) and water density was 


998.5 kg/m?. The data is shown in Figure B.1 


B.3 Mass Flow Rate Correction 

The inlet water temperature from the cooling water sump varies anywhere from 
15°C to 25°C depending on environmental conditions. To account for these 
temperature variations the following function was used to calculate a correction 


factor for viscosity variation with temperature [Ref. 5]. 


Cf=1.0365 - (1.9644E - 3) Tin +(5.2500e -6) Tin? (B.10) 
where: 
Cf = mass flow rate correction factor 
Tin - inlet temperature (celsius) 
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Figure B.1 Horizontal Tube Coolant Flowmeter Calibration Chart 
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The value of Cf for the flow meter calibration (Tin=17.5°C) was 1.0037. 
Therefore, the actual mass flow rate was calculated using the following equation: 


€ 
ALME ES (B.11) 


where: 
M,a = corrected mass flow rate (grams per second) 
ma,7 computed mass flow rate (eq. B.9) (grams per second) 
B.4 Pressure Transducer Calibration 
Three methods of pressure measurement were available on the apparatus: 


1. Direct pressure reading off the Heise solid front - CM-104119 pressure gauge, 
(range 0-15 psia). 


2. Converted voltage readings from the Setra, model 204, Ser. no. 63982 
pressure transducer (range 0-14.7 psia; 0-5 volts; 5V~O psia). 


3. Steam saturation temperature measurement with the apparatus producing 
steam at steady state. The steam saturation temperature/pressure relation 
was utilized via standard steam tables. 

The pressure transducer was calibrated versus the vapor temperature probe 


reading on 12 December 1991. Equation B.12 gives the desired relationship. The 


data is shown in Figure B2. 
P = -2.9360V + 14.7827 (B.12) 


where: 


pressure (psia) 
pressure transducer voltage reading (volts) 


« "9 


B.5 FRICTION TEMPERATURE CORRECTION 
As coolant flows through the tube there is a bulk temperature rise due to 


frictional heating, which is highly dependent on fluid velocity. Although small, this 
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Figure B.2 Pressure Transducer Calibration Chart 
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Voltage (volts) 


Transducer 


can have a significant effect on the calculated overall heat-transfer coefficient. 
Measurements were made for no insert, Heatex, wire wrap and twisted tape inserts 
as shown in Figure B.3, on 5 December 1991. Each data set was curve fitted to a 
third order polynomial which is depicted in Table B.2. Each respective polynomial 
corrects the temperature rise measurement for the heating due to the particular type 


of insert used. 


Table B.2 FRICTION TEMPERATURE RISE POLYNOMIALS 
Insert Type Polynomial 

None Trise = -1.960x10°V* +9.349x10“V? + 1.749x10“*V-2.728x10" 
Wire Wrap Trise — 8.160x10^V? 4- 1.4512x10? V? 4- 2.745x10? V-3.991x10* 


Heatex Trise — 8.160x10^V? - 1.080x10?V* 4- 1.232x10^ V +8.570x10° 
Twisted Tape Trise — 4.070x10^V?-- 4.451x10* V^ - 1.711x10? V-6.440x10^ 


where: Trise — temperature rise (K) 
V = fluid veloci 





84 


(m/s) 


Flow Velocity 


i 
Ф 
v) 
= 
Ф 
Q 
e 
= 
79 
5 
E 
= 
= 
=> 


E No Insert 
B )Heatex Insert 
А Wire Wrap Insert 





со N е LO + со N == © 
© © © = о © © © © 
© © © © © © © © © 


(M) Ө5:їН әлпу}руәӘйшәу 


Figure B.3 Friction Temperature Rise Curves for Heatex Insert, Wire Wrap 
Insert, Twisted Tape Insert, and No Insert. 
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APPENDIX C. SYSTEM INTEGRITY / LEAK TESTING 

As mentioned in section III C, the material for the auxiliary condenser 
penetration plates was changed from aluminum to stainless steel on 24 January 1992. 
The stainless steel screw thread connectors in the aluminum side plate had loosened 
to the extent that leakage could be detected. The cause of this loosening was due 
to thermal cycling of the apparatus and the differential contraction/expansion of the 
aluminum/stainless steel combination. 

An initial leak test was conducted from 20 December 1991 through 2 January 
1992; the results are shown in Figure C.1. The initial mean leak rate was 3.4 mmHg 
per day. 

Subsequent to the structural modification noted above, another leak test was 
conducted 6-19 February 1992; the results are shown in Figure C.2. The mean leak 
rate was found to be 1.7 mmHg per day, an noticeable improvement. In general, a 


leak rate of 2 mmHg per day is considered acceptable. 
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Figure C.1 Apparatus Leak Test I 
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Figure C2 Apparatus Leak Test II 
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APPENDIX D. COMPARISON OF FIXED C; vs FLOATING C; SOLUTION 
METHODS FOR MODIFIED WILSON PLOT DATA REPROCESSING. 

As related in section VI D the "fixed" C; method and "floating" C; method were 
evaluated against the original instrumented tube results in order to choose the most 
accurate method of reprocessing non-instrumented data. The results of the fixed 
coefficient versus modified Wilson plot floating coefficient method comparison are 
shown in Table D. The fixed coefficient method was determined to be more accurate 
than the floating coefficient method 75% of the time. For the three coolant flow 
rates considered for each run, the overall fixed coefficient method mean error was 
+2.0%, and the floating coefficient method mean error was +5.4%. 

For Re <20,000 the error was noticeably higher, which tends to support the 
assertion of a number of researchers that data with Reynolds numbers below 20,000 


should not be used, since the flow may not be fully turbulent. 
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Table D. COMPARISON OF FIXED C; vs FLOATING C; REPROCESSING 


METHODS 


File Name 


FIMAVSHI 


FIMAVSH2 


FIMAVSH3 


FIMAVSH4 


FIMAVSH5 


Flow 
Vel. 


№ 
(m/s) 


4.35 
2.21 
1.16 


4.34 
2.22 
1.16 


4.35 
223 
1.16 


4.35 
223 
1.16 


4.35 
2,23 
1.16 


Film 
AT 


Т, 
(К) 
29.33 


22.81 
18.27 


36.09 
30.39 
24.62 


46.10 
3935 
32.56 


55.20 
47.31 
39.45 


32.83 
28.05 
22:39 


Inst 
Tube 
Data 


h, 
(kW/m^K) 


12.13 
13.18 
13.68 


11.10 
12.03 
12.66 


10.31 
10.89 
11.42 


9.263 
9.869 
10.41 


8.618 
9.589 
10.23 


Fixed 


Method 


11.93 
12.82 
13.16 


11.14 
11.87 
12.50 


10.26 
10.83 
11.40 


9.146 
9.756 
10.39 


8.478 
9.341 
9.865 


To 
еггог 


-1.7 
-2.6 
-3.8 
202 


0.4 

-1.3 
-13 
-0.7* 


-0.5 
-0.6 
-0.2 
-0.4% 


-1.3 
-1.1 
-0.2 
-0.9* 


-1.6 
-2.6 
-3.6 
-2.6* 


Float 


Method 


13.19 
15.17 
17.20 


11.07 
11.72 
12.26 


10.62 
11.44 
12.43 


9.584 
10.52 
11.71 


9.018 


10.36. 


11.66 


90 
еггог 


8.7 
15.1 
25:9 

+16.5 


0.3 
-2.6 
-3.1 
-1.8 


3.0 

5.1 

8.8 
+5.6 


47.5 


File Name 


FIMAVSH6 


FIMAVSH7 


FIMASW3 


FIMASW4 


ҒІМА5УУ5 


Flow 
Vel. 


ЫН 
(m/s) 


4.35 
223 
1.16 


4.35 
2222 
1.16 


4.35 
2:22. 
1.16 


4.35 
222 
1.16 


4.35 
222 
1.16 


Film 
AT 


39:23 
33.76 
27.98 


49.04 
41.25 
34.45 


55.08 
45.50 
36.67 


52.83 
44.76 
36.16 


54.28 
45.39 
36.69 


Inst 
Tube 
Data 


h, 
(kW/m^K) 


9.073 
9757 
10.18 


8.960 
9.685 
10.17 


9.379 
10.15 
10.73 


9.531 
10.15 
10.67 


9.405 
10.08 
10.61 
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Fixed 


Method 


8.953 
9.447 
10.25 


8.939 
9.513 
10.01 


9.384 
10.11 
11.18 


9.558 
10.22 
11.18 


9.453 
10.14 
11.20 


% 
еггог 


-1.3 

-1.3 

0.6 
-0.7* 


-0.2 
-1.8 
-1.6 
SN 


0.1 

-0.4 

4.2 
41.3" 


0.3 
0.7 
4.8 

+1.9 


0.5 

0.6 

5.6 
+2.2 


Float 


Method 


9.084 
9.667 
10.64 


8.637 
9.010 
9.209 


92723 
9.813 
10.62 


9.528 
10.17 
11.08 


9.367 
9.985 
10.90 


% 
еггог 


0.1 

1.0 

4.5 
+ 1.9 


-3.6 
-7.0 
-9.4 
-6.7 


-1.7 
-3.3 
-1.0 
-2.0 


0.0 

0.2 

4.2 
41.5" 


0.4 
0.9 
2.1 

+1.3* 


File Name 


FIMASN4 


FIMASN6 


Flow 
Vel. 


V, 
(m/s) 
4.35 
2.22 
1.16 


4.35 
2:22 
1.16 


Film 
AT 


T 

(К) 
46.24 
34.27 
24.82 


45.93 
34.40 
25:52 


Inst 
Tube 
Data 


h, 
(kW/m?K) 


9.947 
11.12 
12.21 


9.872 
10.84 
11.81 
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Fixed 


Method 


10.07 
11.34 


9.823 
10.86 
14.35 


Io 


еггог 


182 
2.0 


иби 


Float 
Method 


9.512 
10.24 


9.302 
9.698 
11.77 


% 
еггог 


-44 
-7.9 


-6.2 


5.8 

8.7 

0.3 
1:497 


APPENDIX E. UNCERTAINTY ANALYSIS 


In the measurement of a physical quantity, there will always be a difference 
between the measured value of the quantity and the actual value. The magnitude 
of this difference depends on the accuracy of the measuring device calibration, 
operator experience, environmental effects, etc. Eventhough the error associated 
with a single measurement may be rather small, the error may grow to substantial 
proportions when combined with other measured quantities in a given calculation 
scheme. The best estimate of the difference between a calculated or measured 
quantity and the actual value of the quantity is known as the uncertainty. 

The uncertainty may be estimated by the method of Kline and McClintock 
[Ref. 24]. This states that for a quantity R, which is a function of several measured 
quantities (R=R(x,,x,,x;,....X,)), the uncertainty in R is given by the following 
relation: 


W, = (E.1) 








2 2 2 
шала Ен 
ox, Ox, Ox, 
where: 

Wr = the uncertainty of the desired dependent variable 


X1,X7,X3,--..X, = the measured independent variables 


W,, W,, W,,...,.W, = the uncertainties in the measured variables. 


Georgiadis [Ref. 5] gives a complete description of the uncertainty analysis for 
this experiment. The uncertainty analysis program written by Mitrou [Ref. 25] was 


used to calculate uncertainties. The uncertainty for runs using an insert ranged from 
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+2.2% to +4.6%. The uncertainty for runs without an insert ranged from - 10.946 
to +17.1%. Sample outputs of the uncertainty evaluations are included in this 


appendix. 
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APPENDIX F. INSTRUMENTED TUBE CONSTRUCTION 


The medium diameter instrumented tube of Poole [Ref. 6], with six wall 


thermocouples, was fabricated by welding together a copper tube in three pieces 


after imbedding capillary tubes at mid depth in the tube wall. Teflon thermocouples 


were inserted into the capillary tubes for tube wall temperature measurement. 


An attempt was made during this study to construct instrumented tubes of 


small, medium, and large diameter in the following steps: 


UO EI? 


lake thick-walled copper tube bar-stock and machine to specified inside 
diameter, outside diameter, and length. 


Cur four evenly spaced slots, of p depth and width to accommodate metal 
sheathed thermocouples, from halfway along the tube to the end. 


Solder metal sheathed thermocouples into groves. 

Cur copper strips from another tube to fit into the top of the slot. 
Clamp copper strips to the slots (using jubilee clips) and solder in place. 
Turn off excess copper from strips to original outside diameter. 


Send tubes to plating shop and plate the whole tube to a given thickness of 
plate. 


Return tube to NPS machine shop and machine back down to desired outside 
diameter. 


The process was completed through step 2 of the above procedure and the 


tubes are ready for step 3. Several of the latter steps were attempted using a 


practice tube with monel wire placed in the slots on the tube vice thermocouples; this 


was to ensure the process was safe, and to prevent destruction of assets (small, 


medium, and large diameter tubes machined to specifications, and the associated 
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metal sheathed thermocouples). The following summary documents the steps taken 


and lessons learned. 


1. Completed tube fabrication p step 2 of the fabrication procedure 
(above). A practice tube was used from this point to continue the procedure. 


2. Cutting the copper strips proved difficult, so the monel (similar to the 
stainless steel sheath on the thermocouples) wires were silver soldered after 
being pinned in place. No lead was permitted in the solder since copper will 
not plate to lead, but will plate to silver. 


3. The practice tube was sent to a local contractor for electroplating. (Bay 
Custom Chrome, Marina, Ca.) 


4. Two different plating procedures were used on the tube: 


a. First, a cyanide bath treatment was used to electroplate the tube. 
However, this procedure severely scorched the surface of the tube, 
resulting in the return of the tube to the NPS machine shop for repair. 


b. After repair, the tube was again sent off and subsequently treated with 
an acid bath procedure, which resulted in an acceptable tube surface 
finish. 

5. The original thickness of electroplate was not thick enough for our 
application, so another acid bath treatment was performed to bring the 
surface thickness to the desired level. 


6. The final product was suitable for our intended application with the exception 
of the following items which need to be corrected: 


a. Existing voids (from non-uniform soldering) or flaws in the tube surface 
would not fill in to give a uniform thickness around the tube. 


b. The practice wires, tube end pieces, and tube interior were not 
protected properly during the electroplating process, resulting in 
unacceptable copper deposition on critical areas. 

It is recommended that the instrumented tube efforts be continued with 
particular attention to the following items: 


1. Silver soldering is an acceptable method of fixing the thermocouples in place, 


since the thermocouple melting temperature 1s ~1700°C (silver solder is 
applied at ~1100°C). 
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If the voids left in previous soldering efforts persist, then cut copper strips to 
place over the thermocouples to hold them 1n place. 


Ensure that critical areas on the tube are well protected during the 
electroplating process. 
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APPENDIX G. DRPINST PROGRAM LISTING 


The program DRPINST, which was used to collect all instrumented data, is 


listed in this appendix. 


107 


с 


) 


c 


КЕ 


2 


сә 


qs 


T 


113 


— — 


сә 


wa 


Ga 
со 


--- 


оз 


C3 
ال‎ 


c 
Lu. 


Gu» 


LO 
С 
сә 


оз 
uJ 
EI 
кез 


u. 


В- 
U? 
— 
mt 


Q 
Ca 


В- 
оз 


=» 


Ё-4 


C) 


Uu 
— 


ЕЕ. 
C 


(23 


С. 1 
Q2 
Se 
Q- 
а- 
C 


(^J 
U) 


M 


me 


U- 
C 


u) 
Ix 
сә 


a 


Єз 
t- 
-4 
Z 


ва 


LJ 
CI 


С> 
єз 
Ес 
ra 
UD 
C2 
a. 


UD» 


we 


113 
< 


c) 


Ca 
Lu 
H- 
Ll 
= 
ас 
H- 
Сэ 


—* 


— 


-4 


ul 
=. 
ui 
a 
ke 


(23 


TR 


C-4 


(13 
xr 


С) 


ІШ 
<F 


сә 


N) 


a) 
xt 


G3 


r- 


a 
p 


st 


аз 
х 


co 


сә 
un 
сэ 


> ДЕ 
=) 


Isla bC 


и?» 
сә 


—» 


= 
u] 
ДЕ 
2: 


C4 
U) 
c 


M 
ш 
© 


т-- 


AT VARIOUS PRESSUR Soe nn. 


st 
1) 
C2 


— — 


U) 
IN 
с; 


s 


Ep 
== 
11) 


U) 
LI} 
C2 


т--- 


108 


SPACING 


dE ERU EY 
FIN 


NSERT URS 
ЕОГЖНН 220201 1021: SO POSITIONS 


NNED 
ЕО ERT INSTALLED 


N=NONE 


4 


ESSET 
1@=1 . 8mm 
1S=1.Smm 
20-2.0mm 


НЕБЕН TU THE TUBE 


Q5=2.Smm 
RUN NUMBER 


6&7 


1055! 
1065 ! 


e 
Сә ~- O? сә 
СЧ | — Od 
| Lu ! | 
ul сз и] и} 
-- ао CO Ol (= 
C -- со сә О 
e сә ч xt 64 
Г- . Cd ч +2 
. —- = . Q 
a? 1 СЧ C © 
2 ^ | ” v) 
U) <r ^ 111 
-- -- LA -- +? 
ч | | -- | v 
—- — Lu LJ І LJ 4) 
а) + m LA и] — + 
. Ц} —- Г- Г- оз 
IN сй Г- N — Г- == 
m c- Г- ro Г- Cd 4) 
LU») c . . а? . © 
Ш a) — . -- +? 
O 0) | 7 -- | v) 
е. с E > 5 
OD nm) — == — — 
Г- a -- | -- -- n 
ег | Lil | | © 
іЭ -- 133) СЧ ul us 
о + — ч Ln c- Сс 
C1 UJ ч Cd Cd из “з 2 
а) сч оз —- a» U) 1) = L 
. 0) С} . а? Сә + L > 
Ul =- к n) : ж 9 +2 (Л QN 
< (D ra | Cd < (Л С © 
x- K) LO > > E S ч- со = 
n Г = p C- Г- Г- о £ Ф + 
>- Ge С! | | | | +2 —4 
<- СЭ 3 Г- | uj uJ u] ЫЈ (. Cn p 
âZ C) ~- A ^ -- X) x) Cd 4) с 
СО а: ae ға ^ e- ы 1 ro K) +? 4) -а4 G- 
Со A- C? з а -- u) а) 0) с! о -41 Q- = 
а: = -4 — Uu) +4 Od о? ГЭ из £ = 
Q- H- m e- C ^ K) UJ . . . . 4j Cn TW 
C) t-4 ma ЫЈ — а) C- из Г- Г- 74 етіс 
ze O Тапы Cy eos w о? а» | | | | U ээ (- 0 
+ > Pe) zt ua ~- ۰ (13 n 5 s 5 D O 
т Ww) Ul et E- Со 22 got TD- Cd Cd с! Cd 4) E + 
LJ e- ш! а. $ê ÛJ CJ CJ 13 С О? | | | | D 0 + © 
а. -а сі а >- >> lad H- D- . LJ uj ul uJ 4 12 0-4 
>> зоо &- - X Uul Q “Ul U) 1) a) гэ ын u {- — -e 
-- 2. O нч n0: vem NS Zeichn xt Г- С м] Ёс O CE 3 
= СЕ = O ZI F- Ww WI UD t) 007 sp a) 89) оз -4 O m ©) 
(е C) H- Lid (ДЕ Uu» 260) ч, Мм > “o? З.2” тэ gq) L1) 10 ua Ш) 
ва Еа ЦІ а ШІ ЕРСІ 02 (3 Г- оо 30 W . . . ` — Q no 
оғ- а г” WO O O) ÛJ F` ~ Ш из Ш Ш +e E- C C) C1 с! Cd C nn 
-10.26-42 010 ОЗ Г- Ғ- ғ- Ғ- Ғ- ө- “1200 с Б 2 = Cu -- 
WwW O Ei U FE IEA O -- ~4 а) + ШІП Ш < 
Е DS мк) © C- -- Я -- Я -- » -- д <Р (1 U} u) 
(ІЗ іп Г- 02 “. m- Єў wb Be + хм. 0 o syu € x (ME ш) Cy нч 
.. =. DIDI UD) ME erg ee О BS 017172 1 02119621) c 


== .. .. LE ee u) .. ae 


O O O O mi u- °. O < f= lO C= U) C>: UO © U er is 


22523425 
4 
В 


Q 
Е- С “ІЗ 7 - СІС222С2С(52С262 (91. ЕСУ | CO CI F- СМ H- ГІЕ- ГІ Е- * I9 м ud 
ОЭ І СЕ Е зы 6 нч К © t- 
кз кї кч кз кз кз кз юз هسم‎ TT O d. C) d Cd. а es шог: 
> >> > хе << OU dr Hr Tr CC Өй * 4 8 ew uj 
GY CO C) C) C) Cla d. dQ ÛJ A GJ TC U) TC UI CC ÛJ wo xX INI 11 
OO OOO OOOO OCCO OOK OC Сэ б бєл а со йа Єз а IPS n. tn 


YM <ч LD 3 Г- О О?) © -- (4 FO «r U WO D- ISO -— NN TO ODD O O O = CIF) ЄР Ш CD 
сә со со со со со со со со со со со сә со со со со ву со Сб Со Со Со OO O O — -- -- -- -- -- -- 


— T~ g= — — т = — == — ge —- т--- = --- we —- T- =- === æa == -= --- — -— — — — gem — ж-- -- —— = 





Кек ты ко ЛИШ nn 1 


чин 


у 
, 


= = гч 

Jw p 
ETE 
Lil 


(NOT COMPL 


ин 


рош 
Take data or re-proce 


1 


ин 


ЗУГТ ыг ес 


"CY 
? 


Е 
с 


| 


sm 


К 
! 


HET 
da 


[2DTMNIT 
Оз дам е 


(ER ,Fp,Fw? IN ММ Еэ ОГТ ММ 


$ с 
4) -4 
+? ¬ 
41 = 
аз лт 
s uJ 
z 222) 
ы» = кз cn 
оз = IL H? H} 4- 
U3 = U 0 ч E 
.. = (9 — — H- 
2- =з uJ = P Е сэ 
Жо 4) ы i-4 ar, [ey ж = 
.. — H- ЮМ 2 ET 
EI: 4] A- U) [б] [Л -4 
au C сс سم‎ ` se «С = 
os an [c] DK 5 = Жи! +? I-4 
ЕЛ z uJ El s В. 2% ІНЕ 
E .. .. i тїр i ГЕ 
.. с! -А e- С) ЕЗ C 
Ze 4) 4) СЕ 5 > ^ 25 
2. Е 5. z ЕЗ z q- 
м A т d- Цаг Lid _J 
ms +> С H- H- we EN I- 
EN C Ч «qo e 2. и СЭ 
2. е т) 8. C3 [=] + = к= 
-+ z {= 17 u -4 
= о: a) с Es ae к 2: а. 
LJ cn Ц] 1 -- -4 
C Б- ü 0 2" O = zt 
< zz + E. il 1: - U) -4 
d LJ 0j (1: — UJ 19 Lil ээг Ll Li- 
70 ~- je C&C ЕЯ Ш! Е- = -- 
Lid H- .. Бин =) 12501211 Er 2 
> єз " uJ uJ 0) $} C2 > EL: Z) -- 
qd 212 “Цэн F- E Qü- —« 13 4) Е - © ae 
C3 8 +2 E) C) C) iur >- Uu- E- а. СЭ U) 
0 Ел Cs ee 5. [EET Sr B > i 4- (27:12 -4 
~ +? Z о s (0 = Hi u- == GE OE BO = (o T 141 
SIE 41 < > =. © I-- Ш = a) | uJ uJ ~- uJ s > E D. ZI 
БЭЭ ^ -- 27 5- parcum Ec 25 оз E uJ 24111 
ZU. = ul >l =- [me Spe E. e e C Z uJ U) ul |- о z 
= e» 4 Сә — we te €) Lil i4 z^ 64 0f d. “- سم‎ pd 
DS ше ШЫ 22 Є 4) ا‎ ын — | H- Uu» C) — c4 Lil єч Li. 
m FH- Li а- +> = ЫЈ Е = C) В "ml. 2? 
A: = [| aj SE а: 2” ü il. -- C? 4 O ul a: 
ul ».. ee ғ- а- Ca W H- uJ 111 ЫЈ ч Б-я 111 uJ -4 dz uJ 
к- O) O) єз .. = [e З < ے‎ qd 4 H- u) IO = CI „4 в | | 
аз < Со С — 3% eo» F4 e- 2° - a ea C) Li 2" ч uJ fe ав АС 
G3 72805290 | u) © LD) 4 uJ c0 Uu) QJ) 69 uJ u) {С Оу Ж” ere- ul 
C- z ЕГ] = C- =z mu) NIS) = — = 3 т tr z z z = т 4 = as: 
- F- O Ц. 0 u uJ z Lid H- cn 
2: — ie F- -ı OF F- o oa’ H- ы ТЕЗ ee 007 I-- E- © Е- F- F- Ғ- F- Fe 2 b- Oe سم‎ 
ола Ем © ры оо а м. а рн ul © 
LJ bj Q- +- = uJ QA- C) F- e e e ш!) Q- ul Uul O. m UJ an UJ Q- mea Еее ЕО. Е- р a. U) Jy 
РИМА ІРІП УСС БИС УСЕШЕ У ul zt d- xg Od: 0^ 0; zt 9 — ll z mI вэ 


C9) O) ^4 C) OD —-4 0) 7-4 LJ UJ O. O. O. ^» QJ -ı O) НН ОА на CL) 4 O) ~+ O. *41 O) GO. O. O. OQ. DH CO нк F4 нэ uJ u- 


<> №) (а 0) C4 а) —- << Dc €» [) XQ) — st C- P) OM rj Lb O) - «xr c- € -- F0 (3 006) -- 


ба F3 EO O 03 C) -- 
С С СЇ Сз [Лу Ту КУ ч <Р РШ) ЦУ HIS DD O O1 GI TC) CD eoe e = pP CI CY CI CI CI CI O C3 ORA 
ось a a OCA CI GCI CI CECI CI CT CI CI са со л їл CI CE) 
= =e - тээ т-- m- --— 90 — ж-- .- T-- -- = e- --— r — = = . - --- -- -- —- -- =~ == =e --- wo .. .. т= سس‎ -- жээ т-- 25 —- тэ 


110 


е лш ЕВ 


ӨР: 


OUTPUT 


[= 


Od 


255 


ТНЕМ 


IF ljob-l 
БЕР? 


Od 


т-- 


g? 
xt 
CJ 


files 


performed for data in fi 


E EXISTING DORIS TILE О 


NAME OF 
711 7: 


THE 


р 


E son 


le 


3 
+ 


ТЫП USING 


а. 


Nrun-18 


v 


C4 


— 


LJ 
0) 


xt 
аз 
OJ 


we 


DATA 


NUMSERSOE 


u 
0o 
--4 
e 
I 
рс 
C) rs 
(QU 
p 
ч 
4) .. 
—ч € 
e4 md 
U_ i 
< Ц 
= 
2 
оз а: 
m| Ц] 
оэ H- 
оз 2 
q- uJ 
сә Р? 
D- D- 
С4 C4 


m4 
Са 
uJ 
C4 
a» 
Od 


e 


n 


Cs 


V) 
12 
دچ‎ 
ma 


<F 
Lu) 
IX) 


—— 


ra 


== 


4) 
Bo 


+2 
+. 


„ч 
13 
4) 
= 


“- 
О 


ГЭ 
C 


1 
g? 


--- 


c» 


“-- 


Ca 


өз 
1 
sr 


— 


GUT Per 


unheated 


Q- 
O 


CJ 


CJ 
Li 
xr 


m 


xr 


we 


аз 
un 
<r 


KNOWN 


WHEN 


(“- 
Г- 
<F 


.NOUN 


аз 
Г- 
xr 


arca 


resistence based on cutiside 


Wall 


A. 
LJ 
LJ 
a3 


со 


—-— 


19 


m= 


end fin Effect: 


ro 


O 


K 


ә 


um MR 


= ов 
- n d DOTE ON " 


51000 


-24 


D 


.. 


mm үң ии 


( 
\ 


111 


По+! 209 


s 


Y ни” 
7 


мм 


{ 
5 


diameter, nc s n’ Po Coe 


и 1 
1 


ШЕЛІ» 


— 


ж. 
1 


ü- 


Cd 
11) 
1 


— 


1648 


— 


LJ 


N 


TH 


211 
О 
ту 

4 


U. 
--4 


С) 
C) 
U3 


=- 


uJ 
uJ 
O) 


O) 
сә 
u) 


эрээ - 


Uu 


8» 
ч- 


23 


E: 
ын 
C 
a 


ud 
e 
ын 
иј 
Ge 
C 


03 
гэ 
u) 


u 
“- 


23 


C 
H- 


0 
гэл 
‚ч 


< 
U) 
-4 
U) 
u} 
4: 


Od 
ч 
u) 


en 


со 
<t 
un 


--- 


Od 


` 


„> 


0 


ч 


л 
U) 


—- 


с? 


ч 
un 
u) 


кз 


С) 


—» 


uJ 


D- 
из 
(13 


-~ 


ud 
0) 
со 


U) 
u) 


—»- 


uJ 
ale 
Ес 


сә 


g) 
u3 
U) 


узсе 


C4 
D- 
(12 


me 


uid 
Un 


и! 


un 
Г- 
U3 


--- 


FINNED""" 


CH 
cu 


UMENT 


INSTR 


СЭ 
= 
+ 
U) 
2 


Ес 
=] 
0: 
Ше 


a3 
[~~ 
U3 


ES 


аз 


сә 
-4 


Un 


ps 
Жао 
-4 
0: 
Gg 


y=- 


a) 
U3 


C) 
C) 


mJ 


С) 


> 


Cd 
GJ 


[*. 
Г- 


Ф 2 
1) 
4) 
DE 
4) 


ч 
ІШ 
[5s 


кз 
ка 
C- 


Г-- 
D- 


Li 
uJ 
03 


p 
ü) 
5 


4: 
оз 


U) 
4: 


оз 
Gi 


=) 
а. 
[=e 


о 


LIJ 
a) 
Г- 


—- 


”7- 


и. 


жээ 
uJ 
Sq 
p 


э 


Li 


єз 


(13 
u) 
г-- 


4 


и. 


ЕЯ 


С) 


ыл 


Li 


1) 
аз 
M~- 


0. 
uJ 
uJ 
u) 


== = 


сә 
C1) 


C)» 


p 
_3 
d- 
Е= 
252 
C3 


Г- 
С) 
сэ 


dq 
UD 


U) 
ын 
a 


c» 
Г- 


js 
23 
a- 
H- 


СЭ 


<q- 
оз 


U) 
4: 


оз 
Сә 
Г- 


ne 
=) 
а. 
ges 
E 
C) 


u3 
(12 


—- 


uj 
и) 
03 


ÛJ) 


wo 


u) 


(D 


= 
=) 
aic 
mi 


—1 
0 


^ 


a) 


0) 
єз 
D- 


Li 
Es 


> 


(11 


un 


a) 


--— 


ІНЕ 
uJ 
и) 
u) 


a» 
C4 
a» 


m 
0) 


~ 


t-4 


<F 
rm) 
а> 


LJ 
ul 
(13 


Сә 
sr 
а) 


c 


оз 
Q) 


Q) 
— 
Q 
сэ 


Lu. 


Е 


C) 
22 
иј 


gI 
ч 
a) 


CP) 


ay 
{_ 
+? 


.. 


а? 


Г- 


Li 
> 


uJ 


1) 
1] 
a 


112 


EUTPUT 


— 
= 


БЕБІ. 
OUTPUT 


= 
м 


` 


t 
vi 


203:"АВ АР40 А153 


p 
a) 


5 


u) 
a) 


те = 


c 
e; 


"АС 


TO 20 


| 


(Уш 
гу 


g? 
a 


r- 


LJ 
y 


и.) 
n 


єз 
o» 


C 


Ca 


I-- 
I-4 
Ч. 
zt 


г- 
Cd 
eo 


CJ 
I 
ro 


y=- 


(- 


а: 
LJ 
ies 


uJ 


0? 
2 
оз 


H- 


Cd 


C 


H- 


--- 
-- 
х.” 


1 


wre 


ser 


un 
ол 


Bu 
uj OQ. aj 


DJ 4 0. 
GJ ~- CJ 


aw) U} aD 
O) O) O? 





=з 


u. 


t 


D- 
a) 
o 


жеге ШИ 
U Ww 
CJ 
z uj 
Сә го 
O) O) 
сэ 0? 


C4 
ы 


0) 
+2 
у? 


CJ 


сә 
C 


un 
С) 
CJ 


113 


а» 


ro 


со 
Cd 


O} 
c» 
Cd 


Lt) 


гэ 


«x 
3 
= 


єз 


Cd 


С? 
C 


18) 
Od 
& 
С} 


и» 


ra 
Г- 
C 


en 


3 
H- 


an 


23 
ү 


ч 
D4 
C» 
C 


U) 
СІ 
сэ 
C4 


В- 


uJ 


== 


O? 
CJ 
сә 
OJ 


Cd 
m 
сэ 
Cd 


<F 
sr 
C 
Od 


i- 
из 
сә 
Od 


U) 


a) 
и? 
c» 
C4 


со 


42 
U) 


оз 
и? 
С) 
qu 


20 


а. 
uJ 
ul 
(13 


зас 
С) 
С» 


U) 


СУ 


p 


c. 
Cl 


Cl 
а] 


^ 
т 


+? 
92 
0 
— 


Q- 


U) 
Г- 
Сә 
бИ 


(51 
nj 


` 


Q- 


СЧ 
+ > 
v 
T 


+? 


С- 
D- 


C4 


4) 
Q- 


--- 


єз 


ul 


UD 
U) 


Г 
(12 
С) 
YT 


CI) 
SJ 
С1 


(19 
oO) 
С» 
CJ 


Ч) 


a3 


C4 


uU) 


жээ 


C 


[ez 
a 
> 


Г] 


H- 
22 
-4 


а. 
а. 


л 
fey 


—- 


H- 
2° 
-4 


oe 
а. 
№ 


m 


езе 


ШТ! 


tn 


I-4 
{N 
22 


r- 
2° 
4 
8 
N 


сә 
<F 


-- 


(CJ CE CU 


[€) 
< 


Cl 


u) 
<F 


—- 


C1 


t-4 
C) 


ud 


O) 
<F 


== 


Ci 


-- 
é- 
р 


а. 


С) 
и? 


ya 


CJ 


9/ ин 
76 


С 
aj 
i- 


а. 


1) 
а. 


БЭР 
47 
y) 


1 
1 


= 


C 


.- 
4) 
.. 
0 
ЖЕ 


Mclal 


“-” 


C 


aj 
Q- 


U) 


a) 
(3 


— 


Od 


114 


| 


-4 
QZ 
Q- 


U) 


= 


+2 СЧ 


сә 


Г- 


Cd 


а. 


زرا 
)0 


С) 
C- 


DJ 


TREN 


АМЕ 


U) 
D- 


C4 


О. 


aa 


g) 


Г- 
DJ 


В- 


- 
о: 


Cd 
a» 


e 


Od 


ин 


‚stem 


Зу 


um 


ЕЭ 


uJ 
a) 


a» 
а) 


CJ 


PUN (1-үҮ,9-М37",0К 


—» 


Lid 
= 


Ey 


Сэ 


‚ч 


ч 
ce 


— — 


CJ 


CJ 
Od 
CJ 


o 
OJ 
CJ 


тээ 


CJ 
Od 


Ifm-2 


<F 
=) 
CJ 
Od 


<F 


a) 


—- 


uJ 
p 
Е 


Со 
у 


ч- 
-4 


ц: 


-4 


ro 
n 
CJ 
CJ 


u) 
st 
CJ 
C4 


: 


у 
/ 


Л 
a 


ITysyvG@tEmft 


=F 


Te! 


<F 
Cil 
C 


a) 


CJ 
Cd 


1 


س مسو 


ro 


un) 


I 
C 
H- 


O 
— 


1 
Cd 


<F 


—- 


(1 


— 


aO» 


— 


uid 


TOUT2  DELT! DELT2 


TINS 


Od 
Ь- 


TOU 


" 4146 


USING 


(Quartz? 


(Metal?) 


> 


we 


H- 
nm 
‚я 


0: 


<F 
U) 
CJ 
Cd 


C4 


115 


U) 
وسم‎ 


а: 
uj 
P- 
єз 
az 
a. 


a) 
u) 
OJ 
N 


G3 
~- 
с} 
N 


а. 


иі 
а) 


Г- 
С} 
Od 


C4 
fos 
C 
OJ 


un 


C4 
CJ 


с) 
c- 
СЧ 
Od 


~- 
с? 
Cu 
Od 


сә 
сэ 
N 


ul 
3E 
(== 


гэ 
O) 
CJ 


a. 
ul 
ul 
a 


и? 
o) 
od 
С} 


a 
о? 
CJ 
Od 


— 


CJ 
— 


C 


Uu 
С) 
Гэ 
C 


СО 


m 
Cd 


xt 


ГЭ 
C4 


Li) 


ГЭ 
CJ 


п и 


г4 
H- 
! 
nm 
-- 


Ee 


о 
ы- 


о 
OJ 


сә 
Сі 
Гэ 
CJ 


-4 


C3 
Z- 


иі 


N) 
C4 
ro 
Od 


u) 
CJ 
гэ 
с! 


С) 
га 


ай! 


DTSE" ни 


Tanp: 


O) 
ГЭ 
гэ 
С] 


Cy 
Е- 


0 
Cl 
p 


< 
<F 
mM 
CJ 


HA 


— 


C) 


zen 


ul 


-y 


ui 
alc 
Е: 


ro 


a) 
sr 


С: 


РЕН 


Ean TEE 


CO 


27 


f:) 
C4 


rd 


к- 


r4 


س 


Ee 


Сә 
U) 
ro 
CJ 


гэ 
Q 
[== 


O 
Cl 
I 


19 
ГЭ 
C4 


из 


os 


O 
ЕС 


„+ 
v- 
H- 


м... 


cn 


~ 


43 


гэ 
и) 
n 
ы 


G3 


«x 
Q- 
t-4 


xt 
из 
ro 
CJ 


Lia 


N) 
C4 


O? 
М) 
гэ 
C4 


ra 


nur 


Жо 


Cd 
Ч) 
ra 
CJ 


[< 


т 
LI 
[o 
C4 


a) 
u) 
г] 
C4 


ет“ 


`~ 


CJ 


ed 


C 


-4 
Q. 
ч 


4- 
> 


24 


~- 
гэ 
CJ 


EY 
=? 


.75E-4*Uyu £9, 3EE-4*Uu^2-1. 835E- 


u- 


وسم 


гэ 
а? 
ГЭ 
C 


116 


/М%4,45Е-4»/ш”2%4,07Е-С» 0,75) 


Е вы Ее IBE-EWUJ^5) 


(-5,44Ё-541.71ЕЁ-3» 
QQE-4242,75E 


ПЕЕ Сотто" 


Inn! 


E 


T 
n 


O2 
a) 
I 


C4 


=A 
=/ 


а 
ww 


гэ 


о"“Т2с-( 


— 


T2 
FIND OUTER 


10 


TE 


х 


Inn=2 THEN 
QeMd*Cpu*tT2c-T13? 


TE! 
Tumi=Tum 


LF 
Ї 


2413 
417! 
2418 


” 
4. 


nm) 


ч 
CJ 


МРЕНА 


i 
= 


FIND INNER WALL T 


TO 


ERATE 


<F 
CJ 


СЧ 


a) 
Od 


C4 


Г- 
C4 


CJ 


С) 
N) 
<F 


CJ 


С 
Ir 
х 
CJ 


wes 


+? 
ЕЗ 
С 
I-4 


C 
Ес 


c 
к 
QO: 
Li. 
r3 
n 


xt 
Od 


xr 


гэ 


C4 


f= нии 


2a 
к-з 
QZ 
A- 


и» 
+ 
ч 
C4 


c- 
xr 
x 
Cd 


52 
+? 


Q- 
єз 


Сә 


х 
CJ 


MuwskNMuw 


19 
ын 
C 


Uu» 


H- 


а. 


Єз 
a» 
<F 
C 


117 


В 


“. 7 


ч 
а. 


Г.) 
13 


4j 
СЕ 
ч: 


QD 


st 
re 


a 
-? 


E 


4 
ay 
Q- 
L 


U) 


nt 


Ca 


eA 


е 


eo 
u 
<F 
OJ 


=ч 


mE 


Ра 
„4 


IE 


xr 
C- 
xt 
С 


Г- 
Г- 
<F 
Od 


a» 
M- 
xt 
Cl 


0) 
c- 


OJ 


©) 
a) 
xt 


Cd 


CJ 
(9 9) 


С 


г) 
a) 
<F 
Cd 


tf) 
сз 
sr 
CH 


u 
(12 
xr 
C4 


C- 
aì 
<F 
C1 


QJ 
a» 
«t 
Od 


оз 
ад 
<f 
С! 


С) 
g? 
xt 
С] 


„+ 
au 
— 
к-4 
QZ 
а- 


c 
сә 


[d 


e 
13 
% 
+2 
V) 


(n 
ЕЕ 


a 
u- 


en 


Q- 


Е 


Cs 
Uu 
Od 


М 


C4 


— — سم 
,+ — 
— — — — 


ZEN 


117 


C3 
us 


rd 


+14 
C3 
ЕЭ 


Lid 
aJ 


СЭ 
E 


u) 
M- 
in 
Cid 


= 


ч 
(1) 
LI} 
P 


uU) 
p 


Lun 


71 


in 


сә 
© 


— 


O 


` 


ap 


2 


(1 


Па 


Пэ 
un 
D4 


ГЭ 
a» 
Lu) 
nd 


u) 
оз 
un 
p 


un 


- 


о 
л 
PAS 


12 
x 


` 


790 
ЖС 


(1 
o 
111 
Ci 


ro 
ge 
un 
ET 


22 
В- 
E 
(9 


U) 
С) 
U) 
o4 


ч! 
n) 
uJ 


МІ 
Гэ 
uJ 
[4 


a» 
n 
uU) 
С} 


118 


!*DosMufsH£g/(Uv^2*K £ *Tdc P) 
-.728*(K £^8*9.81«HfgsRho£^2/(Muf *DosTdc£))^. 


о 
“ 


CF mG 
a os. 


б? 
n 
u) 
Od 


а Е 
& мы} 


Но" н и 
Ba ауы — 


Qf1 
Deor) 


n 


<F 
u3 
C1 


v v 


ИКЕ 
, 


2 
= 


(M 


ГА 
ы 


3 


у 
4 


Y 
ғ 


C4 


H- 
ж. 


a- 
A- 
U} 


Ч) 
Od 


v) 


EIS 


NuRe 


Nure 


m 
Fa 


"Nuss Ret 


эж 


u) 
8) 
22 
С 
аз T 


uuu 
А 
= 


ӨО «- 
0 04 Uu. 


a 

3 

2) 
z- 
— 


3 
22 
2, 
i 
+? 
п 
а: 


х 


ий)! 


u) 
C 


C4 


un 
C- 
U3 
OJ 


U) 
С- 
U) 
C] 


5.14) 


Эрэ 
Q 
гэ 
m 


м) 
m 


С} 


H- 
pad 
єз 
а: 
A- 


a? 
Г- 
U) 
Od 


4) 
O 


+ 


U) 


3E 


о 
а- 


--4 


О 
о 


ини 


G3 
аз 
u) 
CJ 


аз 
(1 
Ga 


152 
a) 
U) 
Od 


<t 
a 
u) 
С! 


и] 
A3 
c- 
G3 
C4 


Q=N)?" Oks 


uj 
E 
-- 


-- 


u) 
e 


C 


4 


N 
C- 
с! 


nj 


Lf) 
C4 
Г- 
Г 1 


и] 
uJ 
A) 


1) 
и? 
ies 
C 


eg 


uj 
U) 
= 
и] 


G3 
i 
~- 
ШІ 


ro 
Г- 
Г- 
CJ 


nm 
a» 
Г- 
D 


ul 
ui 
а) 


Гэ 
9 
a 
C4 


ro 
u) 
C4 


u 
+? 
O 


—4 


(1 


C3 
r- 


st 


—4 


-i 


жу 
сЭ 
ea 
UD 
өз 
d- 


aD 
гэ 


C4 


C4 


H- 


ul 


ж. 


a» 
u) 
a) 
C4 


DES 


1 (A 
zm xr 


un file 


Cres 


+2 


сә 
G3 
ra 


ul 
LJ 
e? OA) 


— 4 

=з D- 

ч- © 
| Сэ 

С) м} 


гэ 


ps 
= 
+4 
а- 
а. 


гэ 
52 
гэ 


ar 
О 
rot 
„4 
“Ҡ- 
ju. 
OQ. + 
ac) 
-2 
= wud 
ч: 
єз icm 
-- CD 
= © 
в [№ 


QN 
— 
= 





Ce 
H- 


-4 


ed 


(69 
a» 
ня 


о? 
u 


ro 
M- 
сә 
N) 


(12 
Г- 
сә 
гэ 


С. 
a» 


ro 


с! 
0) 


гэ 


СЭ 
uJ 


Un 


т-- 


K) 


Q) 


we 


m 


оф 
<%. 


me 
тз 
—4 
Lu. 


z- 


Q 


—⸗ 


CJ 


—- 


гэ 


a) 


-.- 


-1. 
= 
C) 


xt 
C 


О) 


uJ 
IE 
[== 
сә 
+ > 


4- 
| ян. 


на 


Г- 
Cd 


ra 


1795548,54,2328550498 


се 
ы М7 в 


I") 


ra 


ro 


m) 
гэ 


ro 


а- 


U) 
ro 


т- = 


ro 


D? 
ro 


Га 


[<] 


ч! 


C) 
Е- 


©) 
z. 
QZ 


Q 
u. 


ір 
ч 


--- 


m 


a 
<F 


صو 


гэ 


—r 


pes 
>: 
Lid 


--, 


= 


y- 


из 


т-- 


го 


ч 
U) 


wo 


ro 


С 
аз 
Q- 
ж 
ul 

! 


ü- 
D- 
из 
Ix) 


L. 


сә 
сә 
c 


CJ 
C 
C] 
Q- 
со 
U) 


гэ 


m 


Сэ 
22 


ul 


[х 
цэ 


no 


— 


ul 
21. 
E- 


v- 


+ > 
Q- 
мэл 


u. 


+ 


u) 
u) 


oe 


m 


u) 


ар 
LN 


st 


C 
гэ 


O? 
u) 


سا سپ 


гэ 


с) 


Ш 
М) 
Ga» 
CJ 


аз 


-- 


“м, 


“ 


со 
g? 


Сә 
г 
гэ 


шэ 
Э 
u) 
сә 


гэ 
гэ 


Сә 


a 


а. 


Cd 
f- 


Ix) 


r 


Ш 
Г- 


гэ 


—y 


M e 
Lid 
=: 
H- 
C 
it 


د چ 


єз 


ЕЛ 


11) 
a» 
U) 
Ч 
оз 
ra 


gi? 
li 
«г 


<F 
(1) 


АВ! 


u- 


t-4 


ci 


-» 
ж“-- 


ud 


сә 
a» 


o- 


ro 
0) 


Гэ 


оз 
ар! 


Г) 


دچ 
Qe‏ 
A‏ 


TJ 


--4 


- 


— 
С? 
С? 


04 
С) 
CJ 
ГЭ 


un 
Gs 
C1 
ra 


сә 


صو 


м. 


Сә 
Ш? 


Od 


єз 
Є) 
C4 
Г- 
Г- 
<F u- 
С на 


(Nn ci 
9- 222 
T uil 


Rec 
C -- 
C1 DJ 
IO O 


— 


ж. 
Lui 
21 
БЕ 


+2 


кз 


єз+ 


st 


(4 
гэ 


C4 
гэ 


Cd 
P 


ro 


ЧЭ 
CJ 
гэ 
CJ 


Ge 
Яс 
cn 
Is 
[3 
C] 


C 
Гэ 


Jg 
C1 
CJ 
гэ 


19 


ra 
Г- 
Od 


к- 


2 


CJ 
m 
Od 
ГЭ 


гэ 
чэ 
CJ 
a) 
m 


LN 
ro 


Гэ 


гэ 


ч 
xt 
С! 
ГЭ 


120 


+> 
“- 
e+ 


ч 


~ 


На 


--4 
u- 


- 


Ж- 


19 
Cd 
ro 


m 
19 
tJ 
ro 


Sia) 


ro 


J 
E: 
a) 
[- 
Сч 
a 


19 
C1 
гэ 


а» 
un 
C4 


С 
u) 
C4 
m 


19 
U) 
C4 
гэ 


a» 
U3 
Od 
М 


CJ 


>> 
“- 
4 


113 


st 
Г- 
CJ 
no 


r59 


гэ 
Г- 
O4 


H- 
> à 
~ 


== 


= 
-- 


c- 
С- 
O4 
no 


no 


сә 
199) 
C4 


-4 


C) 
“- 
14) 


ro 
a) 
C4 
nm 


uJ 
а) 
C4 
гэ 


C) 
= 
(88) 


4 


Liz 


a) 
О) 


nro 


et» 


-ф 
> 
єз 


"D 


--4 


“ 


2. 
C 
Q 


un 
a) 
OJ 
га 


— 


c- 


BH 
i 


сә 


w- 
+ 


H4 


a) 
а? 
г4 
ro 


c? 
n 
N) 


19 


r3 
CJ 
+ 
>> 
і- 
H- 
<F 


сә 
no 


C- 
C» 
С) 
ro 


10 


ur 
ж. 

қ 
л 


Pat 
A: 
= 
ul 
из 


Od 


orn 


>: 


а. 
> 
и) 


Г 4 
ШЕ 
ГЭ 


гэ 
K) 


U) 


nr 
n 


aD 
ro 
ГЭ 


oJ 
C4 
m 
n 


кз 


и) 


1 
O4 
no 
n 


— 
uJ 
Цэ 
H- 


а) 
Od 
гэ 
ra 


2) 
ro 
nr 


i- 
Cd 


sr 
г“) 
m 
г) 


ra 


C) 


— yp 


и) 


С) 
ч 
Г 
roa 


Га 
<F 
ГЭ 
m 


a) 
«t 
ra 
nro 


гэ 
ro 


CJ 
U) 


м 


ua 
ГЭ 
гэ 


st 


n 
n 


© 
Г- 
гэ 
<) 


CJ 


ГЭ 


4 
C) 
u) 


U) 
Г- 
гэ 
ГЭ 


== 
и 
SE 
Е- 


— 


> 
kb- 
-4 


а? 
Г- 
ro 
no 


IN) 


CJ 


Cd 
a» 
m 
ro 


г) 
гэ 


a) 
a) 


гэ 


1) 


na 
D- 
к 


H- 


Л 


H- 


аз 
гэ 
nro 


ОБВЕС 


D- 
аз 
ro 
n 


121 


c» 
Сә 
xt 


ro 


UD 
c» 
xt 


ra 


С} 


х 


ro 


СЕ 2258 2 0755 


С} 
О» 


un 


зэ” 


Ba 
С} 


a 


-- 
U) 
xt 


o 
C 


U? 


y- 


n 


-4 
23 
LJ 


a) 


<È 


r3 


c 


LJ 
HF- 
42 
Габ. 


М. 


-4 


С} 
<r 
"m 


ro 


<F 
Cd 
xt 
I 


<. 


ais 
Е- 


CJ 
چ‎ + 
=“ 


u- 


4 


Сә 
гэ 
<F 
ra 


10 


0) 
ro 


i) 


— 


ra 
D- 
C 


H- 


22 


n 
гу 
xt 


ro 


о 


u 
I» 
<F 


ra 


“- 
ЁО» 


(3 
а: 
О?) 
ro 
ro 


-4 


> 


ey 


— QS 


<- 


CJ 
xr 
<F 
гэ 


11 
< 
ч 
rn 


и» 
<t 
rna 


xt 
113 
<F 


гэ 


23 
а. 


<- 


а: 


= 
LJ 
а: 


[`- 
аз 
xt 


ro 


ro 


U) 
42 
xt 
ТЕЗІ 


o) 
u3 


m 


UD 


ro 
[`- 
CJ 


из 


ro 


Cd 
хэ 
H- 


- 
`~ 


Се] 
C- 
ч 


ro 


гэ 
Г- 


Г- 
Г- 
и» 
CJ 
un 


>: 
[e 
C 


c 
аз 


> 
U? 


0» 
nO 


u- 


ma 


С) 


и} 


oO) 
C- 
<F 


n 


(1) 
xr 
18 


сэ 
и] 


С] 
11 
аз 
со 
aD 
na 


<r 
uU} 
<F 
АЙ! 


кз 


C3 


—» 


Lil 


112 
<F 


Lid 


(72 
х 
ГЭ 


ГТ) 
(1 
xt 
13 


H- 


Сэ 
С} 
O? 
Сә 
C4 
m 


4 


L_2 


ul 


аз 
0) 
xr 
ra 


(4 
Сә 
U? 
m 


эр 


2° 
и} 
me 
u- 


113 
©) 


ra 


с) 
С? 
113 
N) 


2 
a. 
>: 
и} 


o: 


из 
ГЭ 


> 


ne 


а. 


~ 


D^ 
| 
C» 
<F 


из 
"m 


m) 


гэ 
tJ 
u) 
179 


0 
OJ 
up 
ГЭ 


Сә 


зэ” 


C 


CJ 
ro 
U? 
m 


U? 
IX) 
UD 
7 


I 


Ya li^ 


a) 
гэ 
U? 
ГЭ 


>: 
Lui 


--р 


<- 


т--- 


и? 
ro 


122 


C) 


= 


<. 


т 


us 


un 
ra 


THEN 


О? 
Uu) 


n 


ro 


Cd 
u) 


гэ 


= 


sr 
ке 


т 


+? 
Q- 
+1 


4 


g) 


m 


Cd 
гэ 
4 
с) 


F- 


Q- 


[^- 
un 
гэ 


а 


С 
гэ 


Са 
ж 
ч. 


ss 


Q- 
e 


Г- 
¬ 
19 
ro 


K 


n) 
u) 


гэ 


С) 
Uu» 
C1 


Q- 
SE 


U) 
a) 
0? 
n 


CJ 
03 
U) 
гэ 


сі 
0) 


5% 


l-- 


ИГЕ 


ЦЕ 
Lid 
С) 


Сә 
u) 
ro 


Єз 
a 


— 


22 


гэ 
u) 
ro 


Єз 


ga 


وسم 


F- 


Ч) 
u? 
м) 


(1) 
u) 
№) 


Od 
C4 
u) 
(<) 


(13 


U) 
C4 
Ч) 
гэ 


(1) 
C4 
u) 
1x] 


ярь ~ 


гэ 
гэ 


оз 


u) 
n 


C 
H- 
C 
a 


<t 
ГЭ 
U3 
ro 


وسم 
C‏ 
Lid‏ 


N) 
u 
ro 


Сә 
<r 
u) 
m3 


ГЭ 
ul 
(1) 
m) 


“а 
Ч) 
аз 
rs 


Су 
Ln) 
un 
H- 


H- 


С) 
C- 
u) 
“а 


un 


Q) 
Ес 


єз 
сс 
u. 


[| 


U) 
ro 


м] 
03 
u) 
m} 


ro 
G3 
Г- 
ra 


гэ 


[^- 
(<) 


ГО) 
гэ 
Г- 
гэ 


1 
U) 
Uu) 
E 


u) 
Uu) 

о 
C 


5- 
C2 
ар 
ro 


Г- 
ГЭ 


С) 


(13 
un 
H- 


js 
Uu) 
(^3 


АР 


U} 


(2 
үзе 


> 


li 
-4 


а: 
C 
u- 


u) 
ГЭ 
[>= 
ro 


[- 
гэ 
Г- 
ro 


ra 


a) 
ч 
[- 
ға 


=> 
эг 
un 


»» 


V) 


- 


» 


^ 


H- 


— 


Lu. 
Li. 
СЭ 
гэ 
19 


c- 
K) 


n 


и? 


мэ 


M- 
-- 


Г- 
un 

Q 
Q 


м 


~ 


z- 
C 
Q 


<È 
Uu) 
c- 
ra 


сә 


ылай 


U) 
к- 


U} 
U? 
Г- 
гэ 


123 


0) 
U) 
Г- 
ro 


U) 


— 


аэ» 
—— 


~ 


a) 
un» 

0 
Q 


— 


22 
C 
Сэ 


C- 
C- 
гэ 


4 


сә 
a) 
10 
E- 


к= 


Ч? 


E 
E 


m4 


а: 
Сэ 
ШЕ 


uU) 


Г- 
ГЭ 


Г- 
Г- 
Г- 
го 


C- 
оз 
C- 
ға 


> 
а. 
> 
и] 


а- 


Г- 
O4 
a) 
n 


> 
nr 


ж 


C- 


u) 
n) 


C 
4- 
п. 


з-” 


и? 
d 
u) 


Q 


xt 
с) 
ro 


Г- 
U) 
aD 
ға 


C- 
C- 
СІЗ 
М 


а) 
Г- 
0) 
ro 


оз 
Г- 
а 
га 


LJ 
ul 
aa 


со 
a) 
a) 
n 


e» 
4) 


A 
T 


C3 
p 


12 


— 4 


ч 


cy 
d) 
‚+ 
U? 
U) 
4: 


Od 
а? 
аә 
nm 


а- 
uJ 
uJ 
(13 


152) 
аз 
a) 
FP 


ал 
a» 
ro 


19 
a» 
a) 
rm 


Сә 
C- 


u) 
кз 


ЫЈ 
Е- 
— 
-4 
LI 
A- 


U3 
с) 
OG) 
Fm 


— 
— 


4 


а- 


C- 
a) 
a5 
ГО 


H- 
-4 
ШЕ 
а. 


a3 
a) 
гэ 
m 


O? 
a2 
а) 
Mi 


сә 
(13 
g3 


ни 


U) 
Жо 
-4 


U) 


H- 
т^ 
на 


ШЕ 
а. 


а) 
(1) 
ГЭ 


uJ 
оз 


Lid 


Cd 
0) 
(12 
ra 


гэ 
gI 
99) 
гэ 


u) 
a) 
с) 
Кз 


Ts 


E 


-”4 


а) 
оз 
(1) 
"m 


С 
22 


— 


<£- 


Q 
= 


li. 


(13 
оз 
OD 
Гэ 


Cd 


GJ 
© 
аз 
[з 


С] 
c 
o? 
ға 


со 
оз 
Г? 


a3 
c» 
O) 
ГС) 


(n 

2 

22 
a- 


a 
u) 
Г- 
Сэ 


ao 
ro 


а. 


ul 
02 


0) 
© 
оз 
гэ 


u 
= 4 


„4ч 


и) 
СО 
a 


te 


ГЭ 


C- 
C 
oO} 
ға 


— = 


ul 
02 
мед 
оз 


га 
0} 
м) 


124 


10. 


bl. 


12: 


LIST OF REFERENCES 


Incropera, F.P. and DeWitt, D.P., Introduction to Heat Transfer, John Wiley 
and Sons, New York. pp 566-67, 1990. 


Van Petten, T.L., Filmwise Condensation on Low Integral-Fin Tubes of 


Different Diameters, Master's Thesis, Naval Postgraduate School, Monterey, 
California, December 1988. 


Guttendorf, M.B., Further Development of Filmwise Condensation of Steam on 
Horizontal Integral Finned Tubes, Master's Thesis, Naval Postgraduate School, 
Monterey, California, June 1990. 


Rouk, P. Some Considerations of Data Reduction Techniques in Film 
Condensation Heat Transfer | Measurements, Masters Thesis, Naval 
Postgraduate School, Monterey, California, June 1992. 


Georgiadis, I.V., Filmwise Condensation of Steam on Low Integral-Finned 
Tubes, Master’s Thesis, Naval Postgraduate School, Monterey, California, 
September 1984. 


Poole, W.M., Filmwise Condensation of Steam on Externally-finned Horizontal 
Tubes, Master’s Thesis, Naval Postgraduate School, Monterey, California, 
December 1983. 


Yau, K.K., Cooper, J.R., and Rose, J.W., Effect of Fin Spacing on the 
Performance of Horizontal Integral-Fin Condenser Tubes, ASME Journal of 
Heat Transfer, vol. 107, pp. 337-383, 1985. 


Wanniarachchi, A.S., Marto, P.J., and Rose, J.W., Film Condensation of Steam 
on Horizontal Finned Tubes: Effect of Fin Spacing, Journal of Heat Transfer, 
vol. 108, pp. 960-966, November 1986. 


Katz, D.L., Hope, R.E., and Datsko, S.C., Liquid Retention on Finned Tubes, 
Department 0 Engineerin Research, University of Michigan, Ann Arbor, 
Michigan, project M592, 1946 


Marto, P.J., An Evaluation of Film Condensation on Horizontal Integral- Fin 
Tubes, Journal of Heat Transfer, vol. 110, pp. 1287-1305, November 1988 


Nusselt, W., "Die Oberflachen-Kondensation des Wassertampfes, VDI 
Zeitung, vol. 60, рр. 541-546, 569-575, 1916. 


Memory, S.B. Forced Convection Film Condensation on a Horizontal Tube at 


High Vapor Velocity, PHD Thesis, University of London, London, England, 
September 1989. 


125 


13. 


14. 


137 


16. 


17. 


18. 


19. 


20 


21! 


22 


2» 


24. 


25. 


Shekriladze, I.G. and Gomelauri, V.I., Theoretical Study of Laminar Film 
Condensation of Flowing Vapour, International Journal of Heat and Mass 
Transfer, vol. 9, pp. 581-591, 1966. 


Fujii, T., Honda, H., and Oda, K., Condensation of Steam on a Horizontal 
Tube -- the Influences of Oncoming Velocity and Thermal Condition at the Tube 
Wall, Condensation Heat Transfer, The 18th National Heat Transfer 
Conference, San Diego, California, pp. 35-43, August 1979. 


Rose, J.W., Fundamental of Condensation Heat Transfer: Laminar Film 
Condensation, JSME International Journal, Series II, vol. 31, no. 3, pp. 357- 
375, 1988. 


Dittus, F.W. and Boelter, L.M.K., Heat Transfer in Automobile Radiators of the 
Tubular Type, University of California Publications in Engineering, vol. 2, no. 
13, pp. 443-461, 1930. 


Colburn, A.P., A Method of Correlating Forced Convection Heat Transfer Data 
and a Comparison with Fluid Friction, Transactions of AIChE, vol. 29, pp. 174, 
1933. 


Sieder, E.N., and Tate, C.E., Heat Transfer and Pressure Drop of Liquids in 
Tubes, Industrial Engineering Chemistry, vol. 28, pp. 1429, 1936. 


Petukhov, B.S., Heat Transfer and Friction in Turbulent Pipe Flow with Variable 
Physical Properties, Advances in Heat Transfer, vol. 6, pp. 503, 1970. 


Sleicher, C.A. and Rouse, M.W., A Convenient Correlation for Heat Transfer 
to Constant and Variable Property Fluids in Turbulent Pipe Flow, International 
Journal of Heat and Mass Transfer, vol. 18, p. 677, 1975. 


Lorenz, J.J., Yung, D., Panchal, C., and Layton, G., An Assessment of Heat 
Transfer Correlations for Turbulent Pipe Flow of Water at Prandtl Numbers of 
6.0 and 11.6, Argonne National Laboratory, Argonne, Illinois, January, 1981. 


Coumes, J.M., Some Aspects of Film Condensation of Steam on Horizontal 
Finned Tubes, Master's Thesis, Naval Postgraduate School, Monterey, 
California, December 1989. 


Marto, P.J., Looney, D.J., Rose, J.W., and Wanniarachchi, A.S., Evaluation 
of Organic Coatings for the Promotion of Dropwise Condensation of Steam, 
International Journal of Heat and Mass Transfer, vol. 29, no. 8, pp. 1109- 
1117, 1986. 


Kline, S.J., and McClintock, F.A., Describing Uncertainties in Single-Sample 
Experiments, Mechanical Engineering, vol. 74, pp. 3-8, January 1953. - 


Mitrou, E.S., Film Condensation of Steam on Externally Enhanced Horizontal 


Tubes, Master's Thesis, Naval Postgraduate School, Monterey, California, 
March 1986. 


126 


INITIAL DISTRIBUTION LIST No. Copies 


Defense Technical Information Center A 
Cameron Station 
Alexandria, VA 22304-6145 


Library, Code 0142 2 
Naval Postgraduate School 
Monterey, CA 93943-5002 


Department Chairman, Code ME/He 1 
Department of Mechanical Engineering 

Naval Postgraduate School 

Monterey, CA 93943-5004 


Naval Engineering Curricular Officer, Code 34 1 
Department of Mechanical Engineering 

Naval Postgraduate School 

Monterey, CA 93943-5004 


Professor Paul J. Marto, Code ME/Mx 3 
Department of Mechanical Engineering 

Naval Postgraduate School 

Monterey, CA 93943-5004 


Professor Stephen B. Memory, Code ME/Me 1 
Department of Mechanical Engineering 

Naval Postgraduate School 

Monterey, CA 93943-5004 


Professor John W. Rose 1 
Department of Mechanical Engineering 

Queen Mary College, University of London 

London El 4NS, England 


Mr. David Brown 1 
David W. Taylor Naval Ship Research and 

Development Center 

Annapolis, MD 21402 


LT. Keith A. Swensen 3 


1260 Spruance Rd. 
Monterey, CA 93940 


121 

















DUDLEY KNOX LIBP..ny 
NAVAL POSTGRADO сг “HOO! 
MONTEREY CA 93943-5101 








— — 










> x 

* = ас 

ES E 22 
С), “а К ЖАД As bo» * 
Яо AMO 
PUT Pr 
He eas i 
R С Р 

— 


хээ ER 
x 
EN ог P XA: 
* * oR S 
` V TRIER Н dc 5 4 
Я — Ж 

AN. RO SCR A | | 
PUE | | | С 1 





































де, 
| ЭНХ M 
Эра — 
а. Sch An — 
SER m оо ШУ 
—* 
эс DUM. 50 M S DU 
SOME PM imd SUN om A AR 
E Sia Er H er e ҚҚ ps ORA à at M ATA MN Mu 
S 15 s it EA Pat, S а A Lae 4 5 шоог 
M VL 
ned H VAS RUND RES — 3 T A TA Чин oa Sh 
215170) M — UD n Ns ҚЫЗ kA , 
да A BEN d у 2768 
р i Ч 7 C 
NS Де yr КОН ^ ОО 1752, 
ХХ ROT 51587 5 ГА УЗ КОНУН. КШ, е” 00 : 
PORUM ee d is NICH» RM үн Қ EE ВАЙ — 1 
57 SENT, OP Me d T BE — > 904 
ER — A D КАЙ » NM — | 7 
— 7; A л 0 КОРЧУ 1 ЯЛ OPES f Ж = ЙЕ 
comm КН | 
ж [D чу nal a 1 А Н ГЕ. Nu 2905 
D 55 E [NUR 24 Л А. M ТН ни ee хэ № " Ы inn rm ЕЗ; 4 4 Ч d 4 ШТ Ты P Us Nor: 
Poi $ 3 Хон ЗАМАН CUM SE E Es сай фара КОШЕ? — 
ne 2. Xv И ё ^ (AER Fk KM м A ag КЕ T TA qui |, б Чой i қ), А Е -— x jT — AY "UE 
ТАН va URS ۹ » ARN Aa ka ? eb N K TRN 28 #? Ний LP IW 4: ТЕҢ T УЗ ч БАБЫ Р ши A — 
м ES ir Я EN te n t ы Гэ — ХУАР че х. ТЭ ES * ын BE af 13 AD Y 5 ы "Р Ч Р зо D une я d 0 di " tt 118 ЛЫН Е қ; $5 р E СЕ цо | 
М | ae 5 Ve AH qe 1 A i e NA ч * : ПАТ ИЕ во. | S тұ a 4 ^ 5 өз 3 ‘ ri N i NES Р 
EE EET Kir од 48) г j ba 1 7% DUM 8 .. ФА 1 Ы » й Ы 7 ш Li Li қ x k 2 2 > 
эр ШУ h sd 2 han i Н ia $t Л Lk ча 4 МУ М RT 3, ач ig қ » — — тә» Py H “ЛЫГ ын 4 я БА d m. Н ін А Н 
N 8 — Dd RC , X 1 Bu F жұз 5 s us хан E М е 3, (3 5 1 P EE үс ^ (0$ im noe Ж a Ue Е > uN " 3 Й — 2 р у hr ; if y? 
. A ЖЕҢЕ) КАКА ў АЛЫ MOS cS уи KA м ЖЭ. TEX, UM. prd UE ak fh Ка Vcg E | ov i Жа 
ы ХА AN M ү}; wt 9— LN MAS ا‎ н — Zu hi MS Nj Я EE —* арз р Ст | Ч р M: wt 2 7 Ж: d 
e 5 R 9 —8 ү — Ёл кү Eri қ T ds ur Ka 1 Ги А 214 ЯЛЖЭЭ 348 i m ' СТ к A мэ 5 = g LET 4 g " ETE RC L x 0 4 5 T П Й . 
БЫ ЛЖ > Qu Ж Бета ХЭ 1 Ex Um A г s RES 1; ЖА -Ж к EM QV BUS T6 — Р OM Ж X ON ы j ; E 
45402; um MEM) BR a ШОН кке n DAE NT КАЙ ТОО НЭВЭГ ji. DE Ж i xm Ч DUM o za: — EE NU: M 
4 A ЖО 12° eor ris р 4 A ic, . СҮ at Mr F 7 Ar 1. * ч L i XD 3 % "LI 5 SD 
mm S — Б zd Du И i Y m er LN 30 5 E f Mayer OI А | 4 АС 5 مره‎ РЫ * E қ E E 5 5, = r : 5 ne A ..4 Ь ЛЖ, К ед А n А 
хо Хх Ss РЬ И AULAM. 4 —— Шу П РЕА 4 П Dir =. Ж as L Ц E ; A L * 
—— Bene 5 Td ЕТЕУ: SOME he A КТ ж a — ert og с цг Хо мек — Et 
X Du ДУ er A — ХАН ЕХ к х er x T J ” — A 2 + АН Ж до NEL. ' А ai n жые, op EI 
— s gs — 1. b M T m ne } Pw. . а x OMS 
x 9 ا‎ Р * n. ҮЕ DE ү, РУ SR eet ey — ЖҚ» 8 it e Ku E: D 232 * "T nod А \' mo - Dr N k E н СЕТ TL -.” 565 | Ш T "e 
P 2, vad Жж ж Ч n 7 2 р Ч 1 Кр 1 Wi 2. Жр ОГУ Tu — % S к 1 ee nra B 2 р te Я Sd | 
ХАРЫГ 4 P A У ^ £ ыы? MEM — К: SCH vip 5 r ap Vus A % | Z t ; DN ЁЛ Ын 3 тол, 3 э K 27) ITE .. Он" - Qe : a Й ө i FE ir ғ le 3 
х V 3 У ү} г} KAKAI ER ДА PW CAS 27% өлу Y SNR MU . *8*— ЕТ Jr б 41 26 220 A 4, W DET "m ree бы Der D N 
Nd PAURA а ШЕКТЕ en Ra Eos DIR KA m AC OE к TR ve ERA IC TA бе: hg OMS А | 
Р. n — ae „Аў KR ду ЫЛЕ Ж T - Fur g^ СУУН? Tu E Hi E u ra, Ч TR 127 ы нв ri 7 pi р ү р АЕ Т р n : { 
—— Ba ай Ar TEMO ҚАҚ Be Am vi — Мэс NEM Ди x т. — — i Бог P do — E a 
i > я ы à J » 2 2 5 5 à m x "E E fr "PT : Tt es : " П Р mm ... F П [ME] E " "è аа 
DER орыг — M. MEA D CK ENTE EUIS DEA NETS To — R E 
D ЖА 4 ҮҘ oe 1% LUN hs B RA ЖЕЛІ MH К КЫ Ды i Vene Т, Ч ab x ЖЫ; уж, 120 i | 
1 ЖЭ) ы ен 4-8 М і RAM к) IS ore TER i à р г 4 я G 4 (14 P^ 45,5 Big? 4 "d LE Н в}, „иу As ? 2 TE ar “| — aie 5 ы H 5 : | i | 
p PER н À —* ҚАЛА 4 50; o 2 2: 1 n " E 2 —— e } г ри of T ' re. 2 3 
ы 5554 i m ҰЯ К LA А P АЧ. Bes fs 10 + — * 4 ХУУ c 28 up П a 2 S П " "Y u р NIA us ee M. DD. 23 " un 
р E — Ке ре uA Sos c MU — 12 АЕ. hy 2 en a TU ч г № ар? awe n ж 4 КЫҢ ais 2 EL Ти” j Di 20) VENE D dr ' \ , ы 
n eek р Я 4 3 ын NN Қ ЖҚА s 3 Т k 2 үл Р ^ | ^ 15, ез ITO : "n Я E NP ae $5 Jo ES ШЕ K 1 ч И TE E P П 
if ШОО рұқ; bises Jt B RC зоож Oi ЛАА, оор NC ЛАН " qu a al E NE ша штэ 
ч MTS 12, TRA АТАРДА ХАЯХ, А Же р С, d 20 ) 5 2 M AEN E 
ni СУ? и ua dig 5 ti і ыу Pu Те Хы Мал d ا‎ ЕЛ ү ay 1 Е X 4 at А я D qu = A yo " КОЛЕН а ИТ Ы. TE men D ” УОЛ 
di He ИЗА, ЖЕ, , : ГА! рр E > 3, "| р ^ д E 3 : LX 4 E a ЛАГ .. Ч ыг 27 х, 4 m au? "ү С s T А А Ц 
ТШ УГ QUA QUE 12 ХА п ү —* Nem POR og a Ei ad i No EA de ес — ТИ. ло зэл 2 2 ХЭМ, 2 Ж wi 
m ae WOO Ул — — Edi ET: —* se K ү“ AT al D 8, ^ quw Vy E ASK P Ж, ы bert d а * ям. ко у E Mee ue E е. а ыз — Е 
М м к —— t Ы Н SUSAN Pd PE n Ho AE M TN қ XT шах | Ж) \ | к ES r nA a Д) * CE н PT ‘ t ^ Қу у TT 7 е QC 
in 140 A НИ — қын "a s HER EA i құланы Mu ХЕЛЕН ч Agen n € NE о. — 
Қ ЖАС i uy X uh Mm n Ай ЫГ ty ! к юэ e A ER En Kn 35 ү N Қ a A р к, г Гей ТО ЦУГ do Г k тұ m о s m Va Ls 2 33 М DE 3 E LOUP ۹ ЯГ Li Ы = D LN . 
ый Ut хор is Iu y СИ K Ын T "m А ГИ T EE Е 4 ї m NE СО со Dr E HT 5 * E CN 2) 217 ae ) А “ 
ш я & р ХА H yj 32 m ZA "A en қ Р CEN ETIN E FL х ve 1 vw. * — AE $ Н A 
| i mE n RE со EAN A Зил i 22:49 Vi EA p LEBER ЭЭН E * ЗИ СЕЕ m od 
FA 1 ВОН (9 3 Y ди, ul ж A e — Ps RES Ж VER 5 ШЕТІ ГЭ Мер ) ну Ч Р Т [ а K ы "T ЦЭ П — .y e e ы : و ھا‎ ne П — ⸗ 
mun = RR КЖ NP La RU en де A T5 EN J м "EN, SII Зэн ТЭД Ж” RT mye Ed с Es. | 
LE 55 " [A KU e UR | к? ғы 8 N а? uL 29,15. 2 PELO Dy 4% PT ^ %) ۴ Жэ UM T ы y TN DET 3 n 7 ^ Ч ЭЕ 2 ши ЫН у 2 E : 
He E Dr MA AO ae НА ЖЫН Mi РА 9 E АХА M REN 3,7 г Қо ИШЕ ү Т 7 A 0 g FEE КО. E б . mE Nus ! Й L^ Ar P 
— Y N, * HHS N ARU ЭРЧ ын UE 13 a T ЛЕГІ "E Are x ا‎ $i 1C ТТ x Т Ч "E m .. . 
қ ыр ҚЫН A х1! Чин TAN m КЫ B "D ; қр; ам ЯГ, и vl 5 eke. ^h ЯГ: A 42%) АШ ‘ = 
Von к Мя MC | i ntis $ m D 2: E ERAT. Un. ЯАГ Ро ! ХЭН СО AD ЙГ dE im UT Б ЖИ ШЕ oe "da пасс БАГ мэ D T А — DE 
ДҮШ МН MU Е M Cun TEE К ШИ ЭБЭ) MR "n i S ELT s ЖА е ыр — Ж бн БЭР Ep ХЭЛ? E E 
ki ent р ЛІК CRUS MM MM SM "Ч GAP E Ар NEN: p d — GO 2 " Rh Ч 203 
— i TAE E КӨЛЕ EE MOT O NOME: d D 0. то 
PO > د‎ is ДЕУІ ің, " » Ina” * ТІН. HP “ "DN ШЕ 3! 27 E. Ч 2. a pu ТЫ — ДО à x ағам 5 M 2 4 
A TUE dU EH и 115% а ns hl А “Их E m A aie Tu м CR A д DOM г; —J и Y Жи 1219 x Eh Ч Т TRU Е н 2 цн 
СОН ҰН МЛ и 4 a“ 4 + 4 34 ИН 647. Нар” А m 2 1 H y [n м à 0 З Үр ы ЯГ 3 He ПН Чэ > М Қаз ex д са ы | US E | H 
B^ A 1 x 9 TA A Yu 44 Lt 7 p n ды: "u P ; Lai "ee (2 $ N n A A A ЗАЛ ь Т ‘ А n 5% p «4 Ё 4 A Ets түз mus " М й P ge чол қ р А u 3 E М 1 D Р миа) К va 
M i a я d TR ну Oi P M 1 йү E 3 MA г as VIA IR Pu " i Гм 4 TE 1 4 Е 1 К П Ға "n П Le эн LRL) — "р Л Ч: iu er ds rh ЖАР, " 
hi: ала» ЗЕ г 45% ! M NEH ML tabi ty р 2 MG J. гү 1 у oa VS A: МЫН 7 E ТЕ АН LO "m : 1 TU ч Et * TENA tt 
\ K US А "ү Ж. "m oo ын сэ қ DEM j АНЕ у? HA — Ў DUE 4 JR "m Er DH Aga А АЫ A — LN E. Ы р п Ш ' 
у 12003 — i Popy a d ANT M n ЭЭ a P^ 3 DA T a à a і 2 AT T wt m 4 И D ад JN Yu oa н tt Й | К и d ER " i 2 '' Е i - LI n — 0 5054 б L n а 2 
АС; NU id 1; pe T make 3* AEN кн n E اوي‎ we Ч! ЖУ : Ч, EE Un og ee PS INIT » р P NN i EN. I ans z 5 
ЖАР» POS) Y 0 N! 2! ул ЙГ n? A BO b: 28 4 104” Qe) 53 1, E — Sal و‎ 4 in 1, E "EH 42 2 u "D — "t iy ы ; de ІНІ”. м Ч НЫ D х, Ж Poe 
d RAE qs к Қалы —— vnb p ЧЕ a, Е کن‎ 17243) o DE nt Тит, та LO non CPU E 21 Pug 
yv д К 51% НИ) hay АЛ А D P E А n Я 4 ЯН 4 A t : М 7 
д ИЧ, КТ 1: HAR г eter NUR لرا‎ e P» B. ч үз; ; ay H Cu HU d ҮҮЛ, y г D d Ше” TP i} 74 A f TED ^ бА | Ж” Е ЖИ | oy ын T ; | 
НС ССТ F — MAE Ф К * уу КЕК 7 Л э ZN ‚mi Ли tet At ЫЛ КШ — 29 С A " —— Яг Цан E 13 or: » .. а А A .. ug» 
ну 44, | ; T LAN mm As m — T Е D ИШ ХҮЧ ' 
* m Эй К n LATE E NK UR 82553 | "АМТ ZEE — e D. R x ont ME шон x , E ДЕ 
n) ы? — ямх 54 Dea ЫН A DEUS "i ЯҢА : кон ЖҮ —* P 107 x — er EE RN: E dn ко. DP" 
pt ке « КЕН HR Ё " ҮЕ: : им ity di ien | Ki : эу үзе РЕА 4 Ж. T LL қ — Ч. № 2 ш, цэг ТЕ Das TAG " B PTS } PES ( er: DOSE 5 
хо 5 37, Е Ад n АДЕ И ark h: е НЫ UN >” ЦӨЛ Өз * Zu р » NUR ы PL LI v П Bi ке r р h Ке | [T] PEE П — Е | їй us Ч р a x j 13 vet ош 
ANA CaN Цэ d — УК * р К W 271 dp. iy M * ЯЛ (793 Het m Ee. ы 18) ud zs 47 нэ TEE Tm E E" 
na nr TA MO ‚йа ы! у ab os ! „Ф УС ЭЭЛ d A EL ° ЧЕ | жә. 195 23 
a НУ — НУ, i Ж E ч А rni А 2) ИА ұжы ED Wut s y UR I — КТІ ЛИЕ. A TN EM ЗӨН 
ri ас м Чон 9 t un — хи He i inh hi 2 — i A SE Ae " “ра C Ж; un Ph ШЕШ сэн УЧ | "n р? : и : 
d А — i ы r vt Й UT “7 224 4 pros on — S 7 Ша —— A = цо 4.. : М 
19 yit nd А M y KU SM ши 7 М Y 1 М na m 00 P. ^ s ri 4 XD D B ds ГА А М" NS " 
m. д % 359 A wie Nr ымы ЫМ mhiy Ав? 4 | * 5 E^ fi "ELT “. } м я Е ET ч Ы t [Ra - t 2 i e A 
a Fons E Ed Yeh? KENT UAR M АЛТ E TIT n ek 14 4 НЕХ 9 Nc REN 2020) EAE F x ХӨ 
— 5 XA Чи, u И» қа } 5 vo! ML yv ҒЫ Ї " “2 * " Ты Н Lm At " , "EL 0 TU i H — E т. П ғ 
9 13) AR ^ — aX uS 128 ГИС: ا‎ EV г, ХЭХЭ цг — ХОГ 
5h: н Жин 3 ДМС " КЫРУ м v mo E nh КЫН Е " Ai "T И bial ОЖ. р LA ЙГ) 7 А 1 DAAE Цул м a М 44 Es ‹ ! + Й 2 ш Й \ г 20 А л. Й 
H к A Aye TD ESEL En: у RE г ВЈ П g ЖЕЛ 2 (ЇЕ 1 % x TD» S П AED te . 
ia NA ан ise eae ын М ЖА KR —* к [А wh, uU к vr. n Чэ. PUE rut 54 "үз .- 4 PME я -a tou | 44: Д. кілт — Ом ое 5 “” 5 ; 
Ho ий НУ He i * Hs Me y Ei өң) M A Ч ШИВ С E T Mey ^ A E Же 4 n е З , ҮРЛЭН Hed E 1177 °+ it ee ae Ы x as 3 ., A 
| MEUS don DEN 8 NE * — v cod AN Ww КО 15 
22 i MN 3 ХОЛ s y 1. НЫ и Ф Ws АЖ — 4% 11 e Е Ы у ! Ч E" d Ч sa es LR j : “|ы Ed у i : (Ж | 3 у йг | М 
M } КИ ца a. ний БА $ Le E ТЕН , а-ы SS : 5 — v x р "o EE E 5 d 2 e "A ИХ: 
Я » jh XY Ч Ус е Mo ы) nm sta 4 И 87 5 "ж ы Н ы е P б ь NE RR d у Ley d 2 
M m m NE ЧЕНИН ШҮҮ ы Be e ы” СМ x. " EUN » Ё i EM ы 2 5) 1222 Mx 
k] r D Ч LL "e E MV , ' w 2 АА, E ї : i 
y А d ti 0 M Ins b И ҒАМ» ЕА L D +! 1 P е 1 1 E . 5 "sr Г 1 
V E * n 2 HE E M ы T te z ТУ; V. БЕ > х di sh ! Д ЯС М i? ї АВАГ Аы — VR А an u ey L Лаг. Й L — 
ЯА! on MCA ae үч ха fn 2 US 4 Ч Я Ж жЕ Bo Я ^ EM шоо — r” ei > ar ve” Е ч bey 3 сс 
das p we yas X ES 19777 — n ^ яр | А Аы асуды no ME 2. — 
1 RH и: E „+? { A m EA ber и} ME AN Ну M ТҰ Y Ч “ N 28 ux 7 & d Р EN g u EL n 3 — " Е E ги m i Y 
nn * , ШАР НЫ CR 15 Н ч мэс ‚› \ р ^ $", t ^! v Ф D Эд i i R Ё СА? р 2 % ШЕ ыг 5 er T А e 2. Pd е в D 5 er 5 de М 
ІК. 7 2 ж Н Ї р K oa 1 à CMM UE А ag 11 ЫР. а . ш 22 E | d are ee Ч : i. | р ‘a E E 2 | 
) 9 A n * wr гч | 3 ) m E ^ Н + E en HU ads РА P. u K Ар EE М | TT "e uH © | ED Хэн | ale A ME i 
ЭНЕН АНИ TAR ХАЛ ік " ишиг ру” год NES A ane” — Ч ы 
А E 4 " Ч F FR A 7 4 5 | е? 4 ۰ [CT] E " 0 Li m » 
Зи 5 ХОТ T. TM ) Eu m E га Эу d 
D 4% un 21 : d r TA Ta da E i k D M | Ч E 0 wee or 0 х i L ы И ки 
E "e : ПАН °* Ж Ы E "n Pi ете Fey к: ——— ы ж. П Мі Ku БОР Эд 
Р D 1: | Эр у [t r^ " ДА. 1 р D E 280 TEE UD ' E | VON б, Vr 
7 — г Т á H 21 Ч 44 TE М DE АЖ j = 2 ^ Ы * i 3 эг . 5 1: й P Tu в 
- D My а | | — ee \ xs E Л р пб еы p Ж s р ТЕ 
TEE w 5 | н : u SENT ы" с 3 [3 DE x FH 10 : 2E 
UT 2 : Е en a A ^ ve * ht ۴ M d 1 n , | =! 
A 4 "ee P T L Bei LI 7 La м : e 
D А A 1 " = "I а П .. U 
Er re. EAT — 
E гі, Қ "n 1 , 
у LI № z ы 2 ; LI 
, т Ы É 
D м 0 Н M 





AG 3 ЧИН ! 
Ял, — UP AP 
(қ D TUNE TG LM 
MONA RS ку "t | 
۶ 
И ' I] us б, HET 
— Үг 
Ж LN М М 
иф 
- 4; 1! 
2 Heu 


x 
M E 
Цэ 
SAO EE Y 
442 (54 M | ^ TE 
542% DS ЁС ams КД ОК n m 
: DU Ж КО ah NUT y 
Ён A os 
n HOM Er EAS 
RR M FN 7 
ЫЗ 1 “ы Н БГ د‎ 
e З ч E — 
y ДОР |! Peta, ҚЫ m 
"t 4 3 E П NE i b á 4 Li 
? AN Ж, | X J $ nud * 
n ха 
Lau 3 
* T A y 
a б 
5, 


OE 
T 
ї 7 4 
— DUE ҮН 
> TON A г E | г” ү МЕ hr. en 
^ ns ri SEN ths EAM 155) AY). ЖЕНЕ 
КЕНІ қ A Jit uu "i MM Н т? EA "T Y Ч) v 
SS a — Ж d > DR № ND UE ў); A : 
AM oh т Л c 2 5 
ч sss 92, EN Мо ae ШШ ih BR NE ; (сұ 
N) к қ р un MA р 
— NAE i ug Бал НЫ А RN х RA ХХН нэ x esr 1 | m n. 
An АЛЖ) NS K — ү) Ж DEA 7 Б КИ d D * ILS 1.4 y лэ 
—8 „Ө * E ж Ha » Ge НЕ 
RU НИЯ ЧАК TEN a T 37 1 ne y 
n AM S ү ЖН P p Е" d СЕДІ ОЧ, У unt үр r4 
SEARS ay ae ag Baan АША DI A 
2 "KI ГЁ CIRC LU M Р! ә 
эм: ues u. Жас RA Lus ШҮ : 22 Бай A н 
EL 3u RA T КЫ — 
> —— «^ 9* 
LM 23 d А, E 
Т Nt 4 KE 
! ya E 
. 


5 КЕМ 
E ыг 




















Ц 
[v 
Г ets et L4 Né 
ix MS а. Er 
АКИ 1 СО 
‚ КО Су. 
“ду Хай ) A) ж 
сойр аа in E PE AER 
i N a A Уран P PRA 
2 af Цэ : MS ia К ЖОАН, 
У 3s D t xe Ч ed A ms Ч 
ER Qr PM P E A MS NOE 
AS Ne i PUE m га ШТ! ft DELE 
1 2244 ШЕ ey X 1, Dr 5 Ци SENA y —2 ЧОНО 4 Шыт.) n 
s Р 4 — 2, ү 2 LAC ый Ч H «ЧАДА Тм M ANY D — HE I e n М, t 
^g feud ES UIN eR ees me ШАО 
$ 2 » 5 ү Unt C HUN Ку i hui 
d RS 5 ын ЧАН: P MUN ns К) К ҒЫ SU AHAS Р р ^ yy! Y" pi^ NEXU d) 
er ! SPEER, ЖД A 4 ШИЙ pi d D EAS EN T ma TOM N 0 
{ DOSE ? — ТЭДН ин ин D MEN БОГ, ЖҰЖ J 1524 I — 
DSL КУР ЯН? Л КН ane 1 i e kd ТАИ ы м a, 
К OU 5. ы (UM m КД? d THER — uy у ib. W Эн uw ET 26 ЦЭН п ... 
А ——X Қ 12 р р м 49 ЖЫ. s T A cms 4 А ОРИЗ Г 
P eia Sn o NN n уа, ҚҰЖ» МЯН u s UTR т Ч со m. S 
М [] Ln H 4 M و‎ Н 34.) 1 Ж СЕТ ME nM PORC Ы . ыр Parr E — — " ЕО 4 
ма: 55 s UN Е ЖУ Ж а TN pet, 2 цасан 
—* с; —— НОМ Kar є un + } ле E HERR P ong il * Ч ЕЕ 23 Чэ uS 2.20 % К П қ ЯМ, А — А / OF Ян ? 
— — ee ЦАР Ч, TE FO UR у! үн — | “сү dm P 4 x жа OR A лы у p 5 ЕЯ П 
eu ыы АА » DOS D WO "m | “y ET oum Ew Фон Ew. EE. 
— n А HY dirus 1. s D 27 FL 423 id D dim t ТЭЭН AL |, 
D — ЕДІ vage s 2 ODE б ы d Jl i HE с? м ЗІ. 
— (d [os И ; * ng зү T E Kan XXL put lage iu — "x ET E 
2254 — s Бо — QUARE TNR. EM ТҮЗЕ E RE nn EL EE b X E ^ | 
он | 
5 = qi — RL i E 9 e ie чо. бый A AMD уу? ^ А 7 — Ar Li E 
Е ХА) АӨЖ» СМ MM D M И Pep HC DA E om ah DL Ж ЖЭ ы р ° 
p я 1; m ШОН —B 10 и НО ғ UM DIES — АС у 7 M T Ж гэ ТЕ Ч — ги 2 5 a E С жы: ^ 5 - TC 
оо ЙГ A V КУЙ A M 4 258 пә, Ну LP F ў 21 ЦАГТ НЬ "d Ы A р % [Н "we 23: d 2 e 22) i 8 Pu T d N Pe P 4 
шин x Qu "E ue 257 т ОЕ an atte rou 200 TE Hn a HN. “2607 
ҚАН 2 0 K —5 гр D — Peu вэ — vf n RS JU Жұ М ДЫ A 27. совин AEE | Pee ax 
(e 201 25 ud TE 2. ЕТЕ ҰҒЫН ТҰСТАС io — о КУЕ я № иас 
» y Я 452 E AT cy 15 5 КУ TW et EF Le М5 ла Жы; I > BEN nr te’ П NS pu " |" ЕЗ TIN "C ) ДЭР 17 4 253) "n ы қ ЖЕШ 5 үл: er А н n Ш ^ ^ 
trib aate аг Дон Mot ие — TEE e XR а 2277 $a FEM — — 
D m T M ҒА VN —* AST HO riy HAM A РЕ Ы 42 L "a b 29 fU | " mà ) 2 1% ЗЫ А ^ "m x N — "x А о ор Б C 
KR 3 Ny Muy злі e 444 "T f." ХУ Lig 2 A و‎ As Қ JUAN Lm — i ^ 5 > М 5 x: ТІ. D " x Р TE БЭЭ " У Д 
ite J УАН n — E E E " D M } A ЯС — UE zy n Айс 
имо RAR an ante АЯЫГ Ы —— d ы * d с^ V Цэг ONES ич Ч, 4! А жар ШЫ * УЯ + А; "m DUE CUM bon 
py EHER " у ХЭС —— MEET de $ реу; н а дог * DEEN € ES P | СУЛ) nes Е | p е i» $ | | | 
we " ас АМ a H 1 Pu IN n — Жы Ж : E М e Р Б ы n “» on 
СУЛ ых — 53 TEAM 52 Yap: te ý Цэ + — . gee) ü [m a 4 — "m A = E — P E . 1 5 Lp zr 
UE ДЫН d — НУ d Н Тр u P Y NS p L ТАГ Ғы ТЕН Y Ж. H P. Ч x ы Н р AE D i e Е Й 
N EN d. УЛУК, $ n i ТҮ ^ й M A m ei ий Ч " Es " БО е К ёс Zum A sn T — Яс zn s в e " 2 55: 
KI EE 4а Mu OLET LA "wes | a WE иж 0 E DERE мр n — КҮЗІ. "IDE DN aoe TEC 
314 P eT yt ҮЭ к. 2 2 П 4 210 КЕ З Г "it 1 Ar > T “oy 2 - нг E A X Ар и П TED x “ А 
n LL] s 9 ял Dd Т гч, у Мэй 3 р MN TI 56 A "M. "m Ы “ ا‎ " Ы 0 яг 4 20" 
к ц aan, in d ! a Adi Цаа 1 р E 2 d E Qu 2 
5% Ж? — m в. ЯЛ IR өр : UD E ман R Е AE ut EE А .. 
Ги Шым CR ae] "m Ed е А са 
ы Он d M Ұз 26 Ч E ec e 52 еу : Е — ko 
EE AO diu. D RC ps P 8-0: i TEE 
uns M i СО PR 34 ^. Цаг —— ы » oni > T s » .. Ы б 
bh — ЫШ LOU PUDE П atk р re ne — EE ue , м 
. > П А) 4 E М E е А П G 
ГЇ “Уны, С КТ H 2 Л Tm ғы ы Е = % ы ЫГ, эн П 
4 M у қ i à Fu ы Ы П 4 4 Ы ) nn й 4... i й | d j 
— 25 А T aoe TE р Нь oe Л P 0 * D 
Je E „' PR NET A Tk d | АТЫС 
ШЕЕ К Т Ls d » р — DE Г d ‹ Ta hi 
Xm PX E. E х ын Mr 7 ER 8 
1 Й ШЕ: Р é e ғ“ 
ТЫ 4 Ку cs Қ ГЕ у EM 
4% vee сод PR РҮН Я 
ares el "IT 
31 me 


